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ABSTRACT 

Galaxy clusters’ structure, dominated by dark matter, is traced by member galaxies 
in the optical and hot intra-cluster medium (ICM) in X-rays. We compare the radial 
distribution of these components and determine the mass-to-light ratio vs. system 
mass relation. 

We use 14 clusters from the REXCESS sample which is representative of clusters 
detected in X-ray surveys. Photometric observations with the Wide Field Imager on 
the 2.2m MPG/ESO telescope are used to determine the number density profiles of 
the galaxy distribution out to r 2 oo- These are compared to electron density profiles of 
the ICM obtained using XMM-Newton , and dark matter profiles inferred from scaling 
relations and an NFW model. 

While red sequence galaxies trace the total matter profile, the blue galaxy distri¬ 
bution is much shallower. We see a deficit of faint galaxies in the central regions of 
massive and regular clusters, and strong suppression of bright and faint blue galax¬ 
ies in the centres of cool-core clusters, attributable to ram pressure stripping of gas 
from blue galaxies in high density regions of ICM and disruption of faint galaxies due 
to galaxy interactions. We find a mass-to-light ratio vs. mass relation within r 2 oo of 
(3.0 ± 0.4) x 10 2 /iMqLq 1 at 1O 15 M 0 with slope 0.16 ± 0.14, consistent with most 
previous results. 

Key words: galaxies: clusters: general - X-rays: galaxies: clusters - intergalactic 
medium - dark matter 


1 INTRODUCTION 


Galaxy clusters are the most massive gravitationally bound 
and virialised objects in the obse rvable Universe , with 
masses up to a few x 10 15 M@ (e.g. Kohlinger & Schmidt 
(2014) measure the projected mass within 200 kpc of 
RXC J1347.5-1145 - one of the most X-ray luminous 
clusters found - to be in the range (2.19 - 2.47) x 10 15 Mg). 
Superclusters of galaxies can be more massive and bound, 
but are unvirialised, e.g. Chon, Bohringer & Nowak |2013]). 
Galaxy clusters are interesting as their contents are reason¬ 
ably representative of the contents of the universe as a whole 
and they contain a population of coeval galaxies whose ap¬ 
pearance is affected by a whole range of astrophysical pro¬ 
cesses and the interactions between the processes. 


* Based on observations collected at the European Southern Ob¬ 
servatory (La Silla, Chile). 


The dominant baryonic component of galaxy clusters 
is the intracluster medium, (ICM), a hot (> 10' K), X-ray 
emitting plasma which contains most of the baryonic mass 
(e.g. Lin et ak| 2003). The gas becomes X-ray luminous after 
being heated by adiabatic compression and shocks during 


cluster collapse (e.g. Gunn & Gott 1972 Kravtsov & Bor- 


gam 


20121. A fraction of the clusters have a centre where 


the density is high and where entropies can be low enough 
that cooling should take place on the order of the Hubble 
time (Fabian 1994), but feedback mechanisms quench cool¬ 
ing flows l ] Bower et al.|[2006| |Fabian||2012| is a recent re¬ 
view), adding additional energy to the initially gravitation¬ 
ally heated gas and causing it to be more broadly distributed 
than the dark matter potential. Cluster mass profile estim¬ 
ates from X-ray data, including ICM density measurements, 
are limited to the region in which robust temperature meas¬ 
urements can be made. For XMM-Newton and Chandra, this 
is typically < rsoo; Suzaku and ROSAT can reach r 2 oo with 
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substantially lower spatial resolution (Rei prich et al.||2013 l 
(A is the ‘halo overdensity’ and ta refers to the radius of 
the volume in which the mean density is A times the critical 
density of the universe). 

Cluster galaxies can be used to probe the cluster envir¬ 
onment to greater cluster-centric distances than gas and can 
be treated as approximately collisionless particles moving in 
the dark matter potential well of the cluster. The well known 
morphology-density relation - higher fractions of elliptical 
galaxies in high density environments like galaxy clusters, 
compared to low fractions of elliptical galaxies in the lower 
density field environment ( Dressler|1984 1 - is caused by the 
ICM and the presence of other galaxies, but a precise de¬ 
scription of the way the different components involved in¬ 
teract is not yet available. The critical processes are: ram- 
pressure stripping, where weakly bound gas is stripped away 


from galaxies by interaction with the ICM (Gunn & Gott 


19721; strangulation, where galaxies are starved of cool 
gas in their haloes - gravitational heating combined with 
the active galactic nucleus (AGN)/wind/supernova feedback 
mechanisms already mentioned lead to ICM entropy which 
is too high for effective cooling and replenishment of the 
cool gas ( Larson et al.|1980 1; harassment - gravitational in¬ 
teractions between galaxies which increase internal energy 
and lead to morphological change (Farouki fe Shapiro|1981 


Moore et al.|199 6 Moon yLake fc Katz|1998 1; and galactic 

cannibalism, where dynamical friction reduces the velocity 
of satellites relative to the central galaxy below the velocity 
dispersion of the satellite, allowing it to be accreted on to 
the central galaxy (Ostriker & Hausman!:19771. Galaxies at 
low cluster-centric distances are preferentially harassed and 
starved due to the high ICM density and more frequent en¬ 
counters with other cluster galaxies, including cD galaxies 
which are often the brightest cluster galaxy (BCG). Faint 
and dwarf galaxies with the weakest dark matter haloes and 
the weakest hold on their gas reservoirs fare the worst in 
this environment such that the population of cluster galax¬ 
ies seen at 2 > 0.4 which includes a large fraction of star 
forming galaxies changes by z = 0 into a population domin¬ 


ated by galaxies with very low star formation rates (Butcher 
fe Oemler||1984[ |Dressler et al.||1994[ |Oemler et al.||1997[ 


Popesso et al. |2006| |Boselli fe Gavazzi 2006]) . 

It should be noted that an X-ray selected sample like 
the Representative XMM-Newton Cluster Structure Survey 


(REXCESS, Bohringer et al. 20071, used in the present 


study, preferentially represents objects with higher ICM 
densities, so we could expect that any ICM-dominated ef¬ 
fects would be stronger in our results than in other, optically 
selected samples (e.g. Carlberg et al.||1997 1 or partially op¬ 
tically selected samples (e.g. Popesso et al. 2004). They also 
preferentially include clusters with deep potential wells and 
thus more time available for galactic evolution. Bohringer] 
et al. (2004) show that cool-cores in clusters must be pre¬ 


served on very long time-scales and we could expect these 
to cause distinctive features in the population of galaxies. 
Conversely, if we assume clusters which show disturbances 
in their X-ray morphology are relatively recent mergers, we 
would expect their galactic populations to be less evolved 
than in other types of cluster and might also show some 
trace of disturbance in the galaxy distribution. 

I 11 this work we study the relationship between the op¬ 
tical density profiles of the galaxy distribution and the dens¬ 


ity profiles of X-ray emitting gas in X-ray selected galaxy 
clusters. We explore the extent to which galaxies and gas 
trace one-another and the underlying dark matter. We also 
investigate how the red and blue galaxy populations are dis¬ 
tributed. We investigate the total mass to optical light ratio 
of galaxy clusters, and measure how this varies with respect 
to total cluster mass and morphology (presence/absence of 
cool-cores, regular/disturbed ICM). By taking into account 
the morphology of our sample, we show clear differences in 
the distribution of different types of galaxies in the centres 
of clusters which have had relatively stable morphology for 
a long period of time - massive clusters or those which 
have cool-cores - when compared to clusters with signs of 
more radical recent evolution - disturbed and non cool-core 
clusters. 

The paper is structured as follows. In Section the 
sample characteristics, X-ray data and optical data are de¬ 
scribed. In Section E we describe all stages of the analysis 
including object detection, classification, optical data cal¬ 
ibration, red sequence fitting, radial profile generation and 
luminosity measurements, along with the results generated 
at each stage. The results are discussed in Section S|4] The 
conclusions are summarised in Section Sj5] 

Throughout this paper, radial distances are measured 
in units of rsoo- The influence of the cluster may extend 
further than this, so we typically use the region outside 
1.5 r 500 ~ ?' 2 oo as the off-ta rget regio n, (r 200 = 1 .51rsoo for 
a concentration C500 = 3.2, ( |Arnaud, Poi ntecouteau & Pratt] 
2007)). Magnitudes M in a given broad band filter (R, V, B) 


are signified by a subscript suffix, AB magnitudes by ab, 
Johnson magnitudes by j and K-correction by k- We adopt 
a flat cosmology where h = 0.7 (and h = 0.7 /170), Sl m = 0.3, 
Ha = 0.7 and Ho = 100 h kms -1 Mpc -1 . 


2 SAMPLE DESCRIPTION 
2.1 REXCESS sample 

The REXCESS sample has been compiled as a galaxy cluster 
sample, representative of clusters detected by their X-ray lu¬ 
minosity and independent of their morphology. The sample 


selection is described in Bohringer et al. (2007). The clusters 


have redshifts between z = 0.055 and z = 0.183 and lu¬ 
minosities above 0.4 x 10 44 h 7 Q ergs 1 in the (0.1 - 2.4) keV 
band. The r 500 region, where the mean density is 500x 
the critical density, plus a region outside where the back¬ 
ground can be estimated are within the XMM-Newton field- 
of-view (~ 30'). The mass range of the clusters is M 200 = 
(1.36 - 10.8) x 1O 14 M 0 . They represent a relatively homo¬ 
geneous population in X-ray luminosity, Lx. 

The 14 objects comprising the subset of the REXCESS 
sample used in this work are tabulated with their key para¬ 
meters in Table [T] and comprise approximately half of the 
complete REXCESS sample. The objects were selected by 
right ascension for ease of follow-up observation scheduling, 
and were observed first. 


2.2 X-ray data 


The X-ray observations are described in Bohringer et al. 
12007). Each cluster was observed using all three detect- 
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Table 1. Overview of the REXCESS clusters analysed in this paper. CC = cool-core, D = disturbed, M = massive. Abell Name, RA, 
Dec, z and Lx are from table 3 in |Bohringer et~ al. (:2 007[ l. r §oOi CC and D are from table 1 in P ratt et al.| ( |200 9;); M 500 was derived 
from rsoo and z using the fiducial cosmology. Mas sive ob jects are those with M 500 above the median of the entire REXCESS population, 
2.95 X 10 14 Mq. ID is used in Figures in Section ij.3.9.4|to distinguish between the clusters in the sample. 


Object 

Abell Name 

RA 

Dec 

z 

Lx 

10 37 W 

^500 

kpc 

-M500 

10 14 M Q 

D 

CC 

M 

ID 

RXC J0006.0-3443 

A2721 

00 h 06 m 03.0 s 

—34°43'27.0" 

0.1147 

1.875 

1059.3 

4.21 

D 


M 

1 

RXC J0049.4-2931 

S0084 

00 h 49 m 24.0 s 

—29°31'28.0" 

0.1084 

1.503 

807.8 

1.84 




2 

RXC J0345.7-4112 

S0384 

03 h 45 m 45.7 s 

—41°12'27.0" 

0.0603 

0.495 

688.4 

1.04 


CC 


3 

RXC J0547.6-3152 

A3364 

05 h 47 m 38.2 s 

—31°52'31.0" 

0.1483 

4.667 

1133.7 

5.53 



M 

4 

RXC J0605.8-3518 

A3378 

06 h 05 m 52.8 a 

—35°18'02.0" 

0.1392 

4.478 

1045.9 

4.26 


CC 

M 

5 

RXC J0616.8-4748 


06 h 16 m 53.6 s 

—47°48'18.0" 

0.1164 

1.597 

939.2 

2.95 

D 



6 

RXC J0645.4-5413 

A3404 

06 h 45 m 29.3 s 

—54°13'08.0" 

0.1644 

7.139 

1280.0 

8.23 



M 

7 

RXC J0821.8+0112 

A0653 

08 h 21 m 51.7 a 

+01° 12'42.0" 

0.0822 

0.673 

755.9 

1.44 




8 

RXC J2023.0-2056 

S0868 

20 h 23 m 01.6 a 

—20°56'55.0" 

0.0564 

0.411 

739.5 

1.28 

D 



9 

RXC J2048.1-1750 

A2328 

20 h 48 m 10.6 a 

—17°50'38.0" 

0.1475 

3.215 

1078.0 

4.75 

D 


M 

A 

RXC J2129.8-5048 

A3771 

21 h 29 m 51.0 a 

—50°48'04.0" 

0.0796 

0.767 

900.6 

2.42 

D 



B 

RXC J2218.6-3853 

A3856 

22 h 18 m 40.2 a 

—38°53'51.0" 

0.1411 

3.516 

1130.1 

5.40 

D 


M 

C 

RXC J2234.5—3744 

A3888 

22 h 34 m 31.0 a 

—37°44'06.0" 

0.1510 

6.363 

1283.2 

8.06 



M 

D 

RXC J2319.6-7313 

A3992 

23 h 19 m 41.8 a 

—73°13'51.0" 

0.0984 

0.937 

788.7 

1.68 

D 

CC 


E 


ors (M0S1, M0S2 and PN), and the mean final ex¬ 
posure after cleaning was (1.4 ±0.7) x 10 4 s for PN and 
(2.1 ±0.9) x 10 4 s for each of the MOS detectors. All ex¬ 
posures were cleaned from times of high background due to 
Solar flares and the PN data were corrected for out-of-time 
events. The mean fraction of exposure time lost to Solar 
flares was ~ 0.35 for PN and ~ 0.25 for MOS1/2. 

Cluster centres were set by finding the density peak of 
the X-ray image on a scale of 8.2" (corresponding to 2x 
the PN pixel width), and all radial distances, r, were meas¬ 
ured from these centres, rsoo values are from IPratt et al.| 
12009 ) and were found through iteration of the r 500 — Lx 
relation for morphologically relaxed clusters (eq. 1 in|Pratt| 


et al. 2009 Arnaud et al. 2007). Electron density profiles 


are from Croston et al. (2008). They were derived from sur¬ 
face brightness profiles using the non-parametric method of 


Croston et al. (2006) which performs a direct deprojection 


based on the assumption of spherical symmetry and a reg- 
ularisation procedure, and involves a point spread function 
deconvolution, rather than fitting of a pre-determined gas 
density distribution to the surface brightness profile. 


2.3 Optical data 


We used the Wide Field Imager on the MPG/ESO 2.2 m 
Telescope at La Silla, which is well suited to XMM-Newton 
follow up as it has a similar field of view (34' x 33'). Each 
set of optical data cover the cluster and a region outside r 2 oo 
which we use for the background assessment. The nominal 
resolution is 0.238 "/pixel and the detector is a 4 x 2 array 
of 2kilopixel x 4kilopixel CCDs. 

Dithered observations with total exposure times lis¬ 
ted in Table [ 2 ] were taken in B, V and R bands (ESO 
filters B/123, V/89, and Rc/162). The raw frames were 
reduced and co-added using eso/mvm (alambic). These 
were aligned (shifted and rotated), then cropped to exclude 
regions where any band was missing using the iraf task 


wregister. sextractor (version 2.8.6; Bertin & Arnouts 


19961 was run on each image and the seeing measured by 
taking the median full width at half maximum (FWHM) of 
objects in the unsaturated part of the stellar locus in the 


FWHM-Magnitude diagram. Bright but unsaturated isol¬ 
ated point sources, with SEXTRACTOR-measured stellarity 
> 0.965 in all three bands, were confirmed as point sources 
by eye and used to calculate convolution kernels with which 
to degrade each set of images to a common seeing for pho¬ 
tometry measurements using the iraf task psfmatch. The 
final seeing was equal to the worst seeing in each set of three 
images. Star-galaxy separation was performed on the ori¬ 
ginal, non-PSF-matched images. 

Examples of stacked, flat-fielded images where stars 
have been excised are shown in Figures mmm and 
|A4| (online material). 


2.4 Subsamples 


We use subsamples of the 14 objects in our catalogue, based 
on their X-ray parameters as given in Tablejl] We use the 
morphological classifications of Pratt et al. (2009 §2.3). 

‘Massive’ objects are those with JV /500 above the median 
of the entire REXCESS population 2.95 x 10 14 M@. 

‘Cool-core’ objects have central electron density 
h(z)~ 2 n e ,0 > 4 x 1 0 2 cm and have central cooling times 
< 10 9 yr ( Pratt et al. 2009, §2.3.1). 

‘Disturbed’ objects are classified based on their X-ray 
centroid shifts w;; Wi = < L/r 50 o for d; the projected separa¬ 
tion between the X-ray peak and centroid in apertures with 
radii in the range (0.1 - 1) rsoo- If the standard deviation 
(Wi) is above the threshold value 0.01 rsoo, the object is clas¬ 
sified as disturbed. A detailed description of the determin¬ 
ation of this morphological parameter is given in|Bohringer| 
etTaT] (|20lol §2.4). 


3 ANALYSIS 

3.1 Source detection and classification 

Source detection and measurement was carried out on PSF- 


matched images (see Section 5 2.31 using sextractor in 
double image mode, with R as the detection band. Groups of 
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Table 2. Details of optical observations. Seeing was measured from the images, and the observation date refers to the date when the 
observations were started. A r > r500 and /1 r > r200 are the areas of the regions outside rsoo and r200 which can be used for background 
estimation. 


Object 

-^exp 

h 

^exp 

h 

Rex.p 

h 

Rsee 

arcsec 

Vsee 

arcsec 

Rsee 

arcsec 

-^date 

Vd&te 

-^date 

Ar>r 5 00 
arcmin 2 

Ar>r 2 oo 

arcmin 2 

RXC J0006.0-3443 

0.50 

0.50 

0.75 

0.99 

0.96 

1.05 

2008-09-24 

2008-09-24 

2008-09-24 

647.6 

416.7 

RXC J0049.4-2931 

0.50 

0.61 

1.25 

0.99 

0.97 

1.03 

2008-09-23 

2008-09-23 

2008-09-23 

711.4 

561.5 

RXC J0345.7-4112 

0.50 

0.50 

0.75 

0.88 

0.97 

1.03 

2007-11-19 

2007-11 -19 

2007-11-19 

602.5 

268.0 

RXC J0547.6-3152 

0.50 

1.00 

0.50 

0.76 

1.04 

0.95 

2000-01-19 

2000-01-18 

2007-11-14 

454.7 

315.9 

RXC J0605.8-3518 

0.50 

0.50 

1.37 

0.76 

0.78 

0.92 

2007-11-27 

2007 11 -27 

2007-11-15 

730.6 

586.3 

RXC J0616.8-4748 

0.50 

0.50 

0.75 

0.93 

1.02 

1.08 

2007-11-15 

2007-11-15 

2007-11-15 

740.9 

572.4 

RXC J0645.4-5413 

0.50 

0.50 

0.75 

1.06 

1.08 

1.19 

2000-01-03 

2000-01-03 

2000-01-03 

765.0 

606.2 

RXC J0821.8+0112 

0.83 

0.58 

0.62 

0.95 

0.98 

1.12 

2007-11-15 

2007-11 -16 

2007-11-16 

633.1 

414.1 

RXC J2023.0-2056 

0.56 

0.67 

1.00 

0.90 

1.01 

1.19 

2008-07-01 

2008-07-01 

2008-07-01 

462.3 

85.8 

RXC J2048.1-1750 

0.50 

0.50 

0.75 

0.78 

0.90 

0.99 

2008-07-01 

2008-07-01 

2008-07-03 

667.7 

526.6 

RXC J2129.8-5048 

0.50 

0.50 

0.75 

1.21 

1.25 

1.40 

2008-06-30 

2008-06-30 

2008-06-30 

574.0 

249.2 

RXC J2218.6-3853 

0.50 

0.50 

0.87 

1.18 

1.15 

1.04 

2008-09-20 

2008-09-20 

2008-09-20 

744.9 

573.4 

RXC J2234.5-3744 

0.50 

0.50 

0.75 

1.11 

1.22 

1.41 

2008-06-30 

2008-06-30 

2008-06-30 

674.8 

477.6 

RXC J2319.6-7313 

0.44 

0.56 

0.75 

1.07 

1.06 

1.15 

2008-09-21 

2008-09-21 

2008-09-21 

733.2 

562.3 


> 5 pixels with values > 2er above the median filtered back¬ 
ground per pixel were treated as objects. We used auto_mag 
in sextractor for magnitude measurements. 

Masks were placed on stars: 

< 14 in the Guide Star Catalogue 1.2 

t 

• and/or with prominent diffraction or blooming spikes, 

• and/or with prominent secondary or higher order re¬ 
flections. 


• magnitude Mgsc 
(Morrison et ahpOOl 


Initial placement of the masks was performed using a modi¬ 
fied version of automask.sh from theli (Er ben et al.|2 005) 
and then finely positioned and scaled by hand. Examples are 
shown in Figure [l] Masking introduces sharp edges which 
produce spurious detections not present when using un¬ 
masked images, so objects which only appeared in masked 
images were ignored. When a bright object overlaid the edge 
of a mask, that object was masked individually to avoid 
blending it with nearby extended objects. Objects from par¬ 
tially exposed regions at the edges of the images - a result 
of the telescope dithering and frame stacking - were filtered 
out of the catalogues. 

The sextractor star-galaxy classifier is a neural net¬ 
work trained on a sample of simulated point-source and 
non-point-source images to return a stellarity s in the range 
0-1, where 1 is a point-source and 0 is not. We found that 
good star-galaxy separation was achieved when the max¬ 
imum stellarity from all three un-degraded images was con¬ 
sidered, and objects with s > 0.965 in one of the bands were 
considered to be stars. Bright objects which were not definit¬ 
ively classified (s > 0.8, Rab < 21) were checked by eye. An 
example stellarity-magnitude diagram with images of some 
objects is shown in Figure [2] The figure shows good dis¬ 
crimination (a value close to either 1 or 0) of objects which 
are obviously point-like or extended at bright magnitudes 
(-Rab < 20). At dimmer magnitudes, classifications become 
increasingly random. Because of the finite thickness of the 
Galactic disc and a finite upper magnitude on the star pop¬ 
ulation, we expect the number of stars to drop rapidly above 
the magnitude at which we can no longer properly classify 
stars. Objects at higher magnitudes which are unclassified in 
the un-degraded images are always assumed to be galaxies. 



zfits/ pixel 


Figure 1. Example masks produced by hand for bright objects. 
The circles have radius 25 ". 


By applying the degrading filter, the signal to noise ratio of 
these objects is increased - the spatially uncorrelated back¬ 
ground is suppressed, whereas spatially correlated objects 
on the scale of the kernel used for degradation are enhanced 
- so they can be detected. 


3.2 Magnitude calibration 


Magnitudes from sextractor were converted into AB mag¬ 
nitudes by making an atmospheric correction and a zero 
point correction. We fixed the B band zero point B zelo us¬ 
ing Data Release 7 of the AAVSO Photometric All-Sky Sur¬ 
vey (APASS DR7) catalogue, which has good coverage of 
almost all of our fields, and sufficient coverage (>30%) in 
those fields where coverage was incomplete. 

Using the observations and zero points from |Ziparo| 


et al. (20121, we generated a calibrated stellar locus. By 
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Rab /mag 


(a) Stellarity vs. detection band magnitude. 



(b) Examples of the objects contained in the boxes in Figure|2a| 

Figure 2. Point source discrimination using sextractor’s ob¬ 
ject classifier. Any objects with s < 0.965 are treated as galaxies. 
Objects in box 1 are stars, objects in box 2 were checked by eye, 
objects in box 3 (Rab > 22) could not be classified and are as¬ 
sumed to be galaxies, and objects in box 4 are clearly extended. 
Stars Rab < 13 do not appear as they have been masked (note 
that Mqsc 7^ Rab)- These examples come from the images of 
RXC J0006.0-3443. 


Table 3. Zero points for each of the observations. We found B zero 
using comparison to the APASS catalogue, and then fitted the 
two colour offsets kb—v and kv—R to a stellar locus from the 
literature in order to constrain V Z ero and Rzero • The large vari¬ 
ation in the zero points is due the atmospheric absorption on the 
observation nights; since the extinction coefficients eM were not 
independently determined during observations it was impossible 
to disentangle extinction from the zero point variation. 


Object 

Rzero 

mag 



V” zero 

mag 



Rzero 

mag 



RXC J0006.0-3443 

24.37 

± 

0.02 

24.17 

± 

0.02 

24.70 

± 

0.02 

RXC J0049.4-2931 

23.93 

± 

0.03 

23.92 

± 

0.02 

24.52 

± 

0.02 

RXC J0345.7-4112 

24.52 

± 

0.05 

24.31 

± 

0.05 

24.87 

± 

0.05 

RXC J0547.6-3152 

24.27 

± 

0.02 

24.14 

± 

0.02 

24.44 

± 

0.02 

RXC J0605.8-3518 

24.38 

± 

0.02 

24.11 

± 

0.02 

24.61 

± 

0.02 

RXC J0616.8-4748 

24.45 

± 

0.02 

24.20 

± 

0.02 

24.71 

± 

0.02 

RXC J0645.4-5413 

24.41 

± 

0.02 

24.13 

± 

0.02 

24.61 

± 

0.02 

RXC J0821.8+0112 

24.46 

± 

0.03 

24.21 

± 

0.03 

24.73 

± 

0.03 

RXC J2023.0-2056 

24.31 

± 

0.03 

24.10 

± 

0.03 

24.65 

± 

0.03 

RXC J2048.1-1750 

24.34 

± 

0.02 

24.15 

± 

0.02 

24.71 

± 

0.02 

RXC J2129.8-5048 

24.39 

± 

0.02 

24.19 

± 

0.02 

24.76 

± 

0.02 

RXC J2218.6-3853 

24.30 

± 

0.02 

24.09 

± 

0.02 

24.65 

± 

0.02 

RXC J2234.5-3744 

24.35 

± 

0.02 

24.16 

± 

0.02 

24.69 

± 

0.02 

RXC J2319.6-7313 

24.04 

± 

0.03 

23.79 

± 

0.03 

24.11 

± 

0.03 

galactic plane, where all of our targets are) and thaf 
could cross-calibrate with similar observations. 

; we 


Conversion of magnitudes and colour gradients to the 
Johnson system (for comparisons with the literature) was 
carried out using colour terms from the ESO WFI web 
pageQ The conversion between raw magnitudes M raw and 
Johnson magnitudes Mj is described by 


Rr&w \ 

1 + Cr 

—Cr 

0 

( -Rj 

Traw I = 

— Cv 

1 + Cv 

0 

V 

JJraw J 

0 

—cb 

1 T cb 

V 



(i) 


for JVI zero the zero point in band M, cm the colour term, 
Zm the airmass of the observation and eM the extinction 
parameter. This equation can be inverted in order to find 
standard magnitudes given raw magnitudes. The conversion 
between raw magnitudes and AB magnitudes is described 
by 


minimizing the offset of the stellar locus in each observation 
set from this calibrated stellar locus, we found colour off¬ 
sets KB — V — -Bzero Hzero and K \/ _ j; — Hzero Rzero- This 


method is similar to that of High et al. (20091. The resulting 
zero points are given in Table [3] 

No attempt was made to correct for galactic extinction 
whilst constructing the stellar locus diagram as we found 
that the assumption that the galactic dust could be mod¬ 
elled as a thin sheet was incorrect, and that bluer stars 
(brighter and typically further away) tended to be more 
reddened than redder stars (dimmer and typically closer). 
Instead, we assumed that the observed stellar locus was in¬ 
dependent of position on the sky (at least away from the 


/ Rab \ / Braw \ / Zr en \ 

Vab = Vaw - Zyey (2) 

\ Bab ) \ B raw / \ Zb es / 

/ Or \ / -Rzero \ 

+ I Ov ) + I V zero I , 

\ Ob / \ Bzero J 

for Om the AB offset for band M. The values of the para¬ 
meters are given in Table [4j and the AB corrections Om 
were taken from the ESO mag2flux too£] 


1 https://www.eso.org/sci/facilities/lasilla/ 
instruments/wf i/inst/zeropoints .html [Accessed: 2013-07-24] 
“ http://archive.eso.org/mag2flux/ 
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6 Holland et al. 


Table 4. Magnitude conversion parameters. The offsets Om were 
taken from the ESO mag2flux tool at http: //archive. eso. org/ 
mag2flux/ 


Band 

cm 

cm 

Om 

R 

0.0 ± 0 

0.070 ±0.010 

±0.23 

V 

-0.13 ±0.01 

0.11 ± 0.01 

±0.14 

B 

0.25 ±0.02 

0.22 ±0.01 

-0.07 


The results in the rest of this paper are based on AB 
system magnitudes. 


3.3 Luminosity function analysis 

3.3.1 Galaxy catalogue completeness from off-target 
observations 

The number counts of field galaxies arise from the summa¬ 
tion of luminosity functions at many redshifts, which them¬ 
selves are influenced by galaxy evolution. We use the num¬ 
ber counts histogram to estimate catalogue completeness at 
different magnitudes, and to compare the relative over- or 
under-density of galaxies in the off-target regions of our ob¬ 
servations to other regions in the sky. 

Metcalf e compiled a set of number count measure- 
ment^lfrom Busswell et al. (20041; Frith et al. (20031; Met- 


calfe et al. (|1991|); 

Jones et al. (1991J; Metcalfe et al. 

(2001); 

McCracken et al. 

2000); Metcalfe et al. (1995 1996 

2005). 


The dataset has an intrinsic scatter of ~ 0.3 dex, correspond¬ 
ing to a factor of ~ 2. We use a 5th order polynomial spline 
fit ^Metcalfe (R) as an empirical shape for our field number 
counts histogram fitting. The number count surface densit¬ 
ies S for all of the off-target regions in the sample are plotted 
in Figure [3] along with ^Metcalfe- 

We assume that the fall-off at the magnitude limit 
can be described by the logistic function £f a ii 0 ff (M) = 


^- M falloff 

1 -f- e w fa.Uo(f 


. This function was chosen since it goes 


smoothly from ~ 1 to ~ 0 over a characterisable distance, 
but we make no claim that it precisely describes the fall-off. 
The 50% completeness limit is at M faiioff in this model. 

We fit the function E = / g ^Metcalfe £faiioff where f s is a 
normalisation factor, to the off-target number counts histo¬ 
gram for each cluster in our sample (measured in region 
r > 1.5r5oo), and the results are given in Table [5] We 
also include the fall-off magnitude as a K-corrected abso¬ 
lute magnitude at the cluster redshift, -RABabsKfaiioff- The 


K-corrections are made using the data of Poggianti (19971, 


assuming that the galaxies are E-galaxies. 

The observed galaxy counts were all between lx and 
2x the empirically determined density, consistent with our 
galaxy clusters occupying denser regions of the cosmic web, 
and within the expected bounds of the scatter from Met¬ 
calfe’s galaxy counts dataset. In addition, the shapes of the 
off-target number counts in Figure [3] do not closely trace 
the curve expected from the literature, suggesting that the 
galaxy overdensity of the cluster extends beyond rsoo- 


5 http://astro.dur.ac.uk/~mn/pubhtml/counts/counts.html 
Metcalfe, (2010) [Accessed: 13-05-2013] 



Rab /mag 


Figure 3. The off-target number count surface densities S for 
each cluster in the sample are shown in grey, along with the ex¬ 
pected number counts from data compiled by Metcalfe in solid 
black (see Footnote §). 


Table 6. Schechter function fitting results, valid for the 
MPG/ESO 2.2 m Telescope WFI filters. 


Band 

M* 

mag 

a 

Xr 

Rab a bs k 

-22.26 ±0.17 

-1.18 ±0.04 

2.01 

Fab a bs K 

-21.79 ±0.24 

-1.29 ±0.06 

3.43 

Bab abs K 

-20.87 ± 0.28 

-1.22 ±0.09 

3.67 


3.3.2 Initial cluster luminosity function analysis 


For each cluster, we measured the on-target luminosity func¬ 
tion in the region r < rsoo in the R band, subtracting the 
off-target number counts histogram measured from the re¬ 
gion r > 1.5 rsoo- Each function was normalised by the mean 
density in —21 < R < —17. We truncated each function well 
below the 50% completeness limit, at M f a n 0 ff — 41Tf a ii 0 fr as 


fitted in the off-target region (see Section S 3.3.11. We fit¬ 
ted a Schechter function <f>(L) = <f> exp f° r L» 

the Schechter luminosity (with the corresponding magnitude 
M*), a the slope parameter and <f> a normalisation factor to 
the mean luminosity function. The function and its fit for 
the detection band are shown in Figure [4] and the fitting 
results for all three bands are given in Table [6] 

Given the detection band Schechter magnitude ith, we 
define three groups of galaxies in each cluster: ‘bright’ galax¬ 
ies are those satisfying (R < R, + 2.5), i.e. (L/L* > 0.1); 
‘faint’ galaxies satisfy ( R * ± 2.5 < R < R* + 5), i.e. 
(0.1 > L/L , > 0.01); and ‘dwarf’ galaxies satisfy 
(R > R t + 5), i.e. (L/L* < 0.01). Additionally, we define 
the symbols Lh+ 2.5 = R* + 2.5 and 7?.*+ 5 = R» ± 5. 


The dwarf galaxy population (described in, e.g. Popesso 
et al.| 2006) is not always obvious in our data due to different 
image depths in different observations. The detection band 
data are sufficiently deep to reach the faint limit, R * + 5, 
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Table 5. Off-target number counts fit results. l?ABfalloff is the detection band falloff magnitude in the observer’s magnitude system 
and Rab absK falloff is the corresponding K-corrected absolute magnitude. / g is the normalisation factor by which the literature number 
counts function was multiplied. 


Object -RaB falloff ^ABabsK falloff Wfalloff / E fl A B 




mag 



mag 



mag 






RXC 

J0006.0-3443 

24.07 

± 

0.03 

-14.56 

± 

0.03 

0.153 

± 

0.014 

1.49 

± 

0.06 

RXC 

J0049.4-2931 

23.10 

± 

0.03 

-15.40 

± 

0.03 

0.157 

± 

0.014 

1.33 

± 

0.05 

RXC 

J0345.7-4112 

24.40 

± 

0.04 

-12.76 

± 

0.04 

0.150 

± 

0.011 

1.16 

± 

0.06 

RXC 

J0547.6-3152 

23.34 

± 

0.03 

-15.90 

± 

0.03 

0.150 

± 

0.014 

1.39 

± 

0.06 

RXC 

J0605.8-3518 

23.99 

± 

0.03 

-15.10 

± 

0.03 

0.154 

± 

0.011 

1.53 

± 

0.06 

RXC 

J0616.8-4748 

23.74 

± 

0.02 

-14.93 

± 

0.02 

0.150 

± 

0.008 

1.37 

± 

0.05 

RXC 

J0645.4-5413 

23.74 

± 

0.03 

-15.74 

± 

0.03 

0.151 

± 

0.011 

1.17 

± 

0.05 

RXC 

J0821.8+0112 

24.07 

± 

0.03 

-13.79 

± 

0.03 

0.146 

± 

0.012 

1.27 

± 

0.05 

RXC 

J2023.0-2056 

24.34 

± 

0.03 

-12.67 

± 

0.03 

0.158 

± 

0.009 

1.26 

± 

0.05 

RXC 

J2048.1 -1750 

24.25 

± 

0.02 

-14.98 

± 

0.02 

0.142 

± 

0.009 

1.46 

± 

0.05 

RXC 

J2129.8-5048 

24.02 

± 

0.04 

-13.77 

± 

0.04 

0.152 

± 

0.012 

1.10 

± 

0.05 

RXC 

J2218.6-3853 

22.25 

± 

0.03 

-16.87 

± 

0.03 

0.159 

± 

0.014 

1.19 

± 

0.05 

RXC 

J2234.5-3744 

24.11 

± 

0.03 

-15.17 

± 

0.03 

0.150 

± 

0.010 

1.39 

± 

0.05 

RXC 

J2319.6-7313 

21.93 

± 

0.05 

-16.35 

± 

0.05 

0.143 

± 

0.014 

1.71 

± 

0.10 



Rab abs K /mag 


Figure 4. Mean cluster luminosity function in the detection band 
and best fitting Schechter function model. Before stacking, the 
off-target number count density was subtracted, and each profile 
divided by the remaining object count in —21 < Rab abs K < —17. 


so in the rest of this paper we use the ‘bright’ and ‘faint’ 
galaxies, ignoring the dwarfs. 

By imposing a magnitude cut at f?„ + 5, we lose a frac¬ 
tion of the total luminosity in each cluster. By integrating 
the fitted luminosity functions beyond the cut and extrapol¬ 
ating for very faint objects we can estimate the ratio between 
the total luminosity in galaxies we observe and the total in¬ 
tegrated luminosity in the luminosity function. Luminosities 
based on the ‘faint’ cut need to be increased by 2.5 % in R, 
4.3% in V and 3.1% in B. 

A more exhaustive assessment of the luminosity func¬ 
tions, as well as total luminosity measurements informed by 
the count density profiles is given in Section §3.9| 


3.4 Catalogue contamination by misidentified 
stars 

We estimate an upper bound on the number of stars at each 
magnitude in our catalogues by fitting star counts with re¬ 
spect to magnitude at magnitudes where star-galaxy sep¬ 
aration is robust, and extrapolating this to higher mag¬ 
nitudes. (These count estimates are upper bounds since the 
star counts drop off at a faster rate at higher magnitudes.) 
The number of stars which are expected given the power 
law, but not seen, are assumed to have been misidentified 
as galaxies. The contamination fraction K measured in each 
catalogue at selected magnitudes is given in Table [7] 


3.5 Red sequence selection 

The red sequence is a line in colour-magnitude space around 
which elliptical galaxies in clusters tend to scatter (e.g. 
Bower, Lucey & Ellis 1992| Valentinuzzi et al. 20111. The 
origin of the red sequence and its relation to the mass- 
metallicity relation was explored in a seminal paper by EE1 
imoto & Yoshii (1987). The position of the line changes with 


redshift, and can be used to detect new galaxy clusters and 
estimate cluster redshifts (e.g. Gladders fc Yee|20 00). Spiral 
galaxies in clusters tend to be bluer than the red sequence 
and migrate on to it as star formation fades. The scatter of 
field galaxies in colour space is usually much larger than the 
scatter of the red sequence galaxies and in this study we use 
this observation to increase the signal to noise ratio of our 
radial profiles. 

We fit a line to the red sequence described by 
C mo dei (M) = g (M — Z) + k where C is a colour, g is the 
gradient, M is a magnitude, and k is the colour at a ‘pivot 
point’ Z. The distribution of cluster galaxies perpendicular 
to the red sequence comprises a red sequence of galaxies with 
a relatively narrow scatter (< 0.05 mag in this case) centred 
on the line, and a blue cloud of galaxies with a larger scatter 
(~ 0.5 mag) centred some distance below the red sequence. 
We model this as a distribution of the form 


fit ((/residual) — t-V & 


( C reri<] nal ~°r ) 2 _ +rcBidua l_ °b + 

+if>be ^, (3) 


with CVesiduai (C, M) = C— C mo dei (M), where Or is the offset 


© 2015 RAS, MNRAS 000 , 000-000 



























8 Holland et al. 


Table 7. Contamination fraction for each catalogue measured at the bright limit, the faint limit, and the 50% detection limit. 


Object 

Kr<R*. |_ 2.5 

k r<r, +5 


RXC J0006.0-3443 

-0.00 + 0.05 

0.004 + 0.014 

0.003+ 0.005 

RXC J0049.4-2931 

-0.00 + 0.08 

-0.035 + 0.016 

-0.002 + 0.008 

RXC J0345.7-4112 

0.00 + 0.27 

-0.01 ± 0.05 

-0.005 ± 0.004 

RXC J0547.6-3152 

-0.00 + 0.15 

0.112 + 0.027 

0.154+0.017 

RXC J0605.8-3518 

-0.00 + 0.14 

0.268 + 0.013 

0.319+0.004 

RXC J0616.8-4748 

-0.00 + 0.21 

0.208 + 0.022 

0.270+ 0.007 

RXC J0645.4-5413 

0.02 ± 0.09 

0.331 ± 0.010 

0.365 ± 0.007 

RXC J0821.8+0112 

0.00 + 0.18 

0.223 + 0.025 

0.324+ 0.004 

RXC J2023.0-2056 

0.0 + 0.5 

0.15 ± 0.09 

0.474 ± 0.004 

RXC J2048.1 1750 

0.13 + 0.05 

0.503 + 0.005 

0.5358 ± 0.0033 

RXC J2129.8-5048 

0.00 ± 0.32 

0.05 ± 0.06 

0.183+ 0.005 

RXC J2218.6-3853 

-0.00 + 0.15 

0.368 + 0.011 

0.277+ 0.021 

RXC J2234.5—3744 

-0.02 + 0.04 

-0.057 + 0.011 

0.017+0.006 

RXC J2319.6-7313 

-0.00 + 0.20 

0.06 ± 0.04 

0.222 + 0.015 

Mean 

0.01 ± 0.06 

0.155 + 0.010 

0.2241 ± 0.0025 


of the distribution of red sequence galaxies from the line 
C mo dei, o> the width of the red sequence, and ip r the density 
of the red sequence; Ob, Ob and ipb are equivalent quantities 
for the blue cloud. For well fitted red sequences, o r / a r should 
be small (<C 1). We also define 


w rs 


^residual O r 


(4) 


which is the scaled displacement of a point from the red 
sequence in units of the width of the red sequence; galaxies 
with —3 < w IB < 3 are taken to be on the red sequence and 
those with w IS < —3 are blue cloud objects. 

We consider red sequences in the (R, B — V) = ( M , C ) 
colour-magnitude space (in the AB system), and fit to the 
galaxies satisfying R < 20 and r < rsoo. 

The fitting procedure was: 


(i) Fit of the straight line C mo dei to the colour-magnitude 
points, using a least squares method. (On the first iteration, 
this step was skipped, and we assumed a gradient of —0.044.) 

(ii) Compute the histogram of C re aiduai for all points 
within the magnitude limits. 

(iii) Fit of the scatter model T to the histogram. 

(iv) Select points |w rs | < 3 and use these as the input for 
the straight line fit |(i)| 

Steps |(I)| to |( iv) | were repeated enough times that the results 
stabilized (~ 5 iterations was usually sufficient). In cases 
where the solutions oscillated, the solution with the narrow¬ 
est ay was selected. The fit results are shown in Figure [5] 
To compare our red sequence parameters to those in the 
literature we refitted g and k in (R, B — R) and (R, V — R) 
using the same method, before using 0 and 0 to trans¬ 
form g and k into the Johnson-Cousins (Vj, B, j — Vj) colour- 
magnitude space. K-corrections were made assuming that 
red sequence objects are elliptical galaxies, using values from 
parameters are shown in Table [8] 
11, fig. 3) give values for the K- 
corrected red sequence normalisation fcjK and gradient gj 
for 72 clusters in the (Vj, Bj — Vj) colour-magnitude space. 
The mean values are fejK = 0.87 ± 0.06 (at I4bs = —20) 
and Tjj = —0.044 ± 0.009. These values are in reasonable 
agreement with our results in Table [5] 

The red sequence for RXC J2319.6-7313 is unclear and 


Poggianti ( 1997||. T I 


Valentinuzzi et al. I 


achieving a reasonable fit for this object depended on the 
selection of the solution with the narrowest oy. This is the 
target with the most southerly declination in the REXCESS 
sample, is only visible from ESO La Silla at relatively high 
airmass (~ 1.5) and the observation image depth was some¬ 
what lower than for most of the other targets (as shown in 
Table [5|. We suspect that the uncertainties on the measured 
colours may be underestimated, making the red sequence 
more difficult to see. Before finding the fitting method de¬ 
scribed above, we investigated several alternative fitting al¬ 
gorithms in an attempt to get robust fits for RXC J2319.6- 
7313. Of particular interest was fitting in several bands sim¬ 
ultaneously, i.e. iteration of line fitting in (R, B — R, V — R) 
space with an iterative cut based on perpendicular distance 
from the line. We found that there is a substantial popula¬ 
tion of galaxies which are identified as ‘on the red sequence’ 
using the single band-pair method, but which are excluded 
from it in the multi-band method. Apart from scatter in 
the colour measurements, possible reasons for the offsets 
in the other band pairs are that the galaxies may lie at 
a substantially different redshift, have substantial dust at¬ 
tenuation or AGN activity, or be influenced by the spread 
in star formation histories. The improved discrimination of 
cluster and background galaxies in this case would lead to 
lower background levels overall and higher signal-to-noise 
ratios on the radial profiles. However, this method required 
excellent starting values, and was particularly susceptible to 
blue objects at fainter magnitudes. For RXC J2319.6-7313, 
with its broad red sequence which is poorly separated from 
the blue cloud, the method failed. 

There is no significant change in the the red sequence 
results here - or the subsequent results based on the red 
sequence selection - for the other clusters when using the 
other fitting algorithms we investigated. 


3.6 Radial density model 


We assume a generalised NFW model for the radial density 
profiles, namely 


p(r) = 


Po 




(5) 
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Figure 5. Fitted red sequences for all the clusters in the sample. The red sequence best-fitting line is shown as a solid line, and the 
dashed lines show ±3ru r s- The vertical lines show the Schechter magnitude, R*, and the bright and faint magnitude limits, R*+ 2.5 and 
R „+5 respectively. The star denotes the BCG. Objects below the lower dashed line are ‘blue,’ and objects between the dashed lines are 
‘on the red sequence.’ 


Table 8. Red sequence parameters in the Johnson-Cousins colour-magnitude space (Vj. Bj — Vj), where gj is the gradient, and fcj = 
Vj, .Bj — Vj at V a bs = — 20. Icjk are compensated for evolution and redshift dependence on colour using values from Poggianti (1997) 
assuming that the red-sequence comprises of E-galaxies. 


Object 

93 

fcj 

kjK 



mag 

mag 

RXC J0006.0-3443 

-0.063 + 0.024 

1.12 ± 0.16 

0.74 + 0.16 

RXC J0049.4-2931 

-0.035 ± 0.031 

1.23 + 0.16 

0.87 + 0.16 

RXC J0345.7-4112 

-0.053 + 0.068 

0.995 + 0.155 

0.84 + 0.16 

RXC J0547.6-3152 

-0.059 + 0.015 

1.18 + 0.16 

0.71 ± 0.16 

RXC J0605.8-3518 

-0.077 + 0.017 

1.11 ± 0.16 

0.66 + 0.16 

RXC J0616.8-4748 

-0.045 ± 0.022 

0.963 + 0.150 

0.57 + 0.15 

RXC J0645.4-5413 

-0.071 ± 0.012 

1.15 ± 0.16 

0.64 + 0.16 

RXC J0821.8+0112 

-0.028 + 0.010 

1.11 ± 0.15 

0.85 ± 0.15 

RXC J2023.0-2056 

-0.038 + 0.042 

1.10 + 0.15 

0.92 ± 0.15 

RXC J2048.1-1750 

-0.064 + 0.021 

1.03 + 0.16 

0.56 + 0.16 

RXC J2129.8-5048 

-0.032 ± 0.009 

1.12 ± 0.15 

0.86 + 0.15 

RXC J2218.6-3853 

-0.055 ± 0.059 

1.22 ± 0.19 

0.77 + 0.19 

RXC J2234.5-3744 

-0.066 + 0.007 

0.997 + 0.143 

0.51 ± 0.14 

RXC J2319.6-7313 

+0.007 ± 0.034 

1.15 ± 0.16 

0.83 + 0.16 

Mean 

-0.049 ± 0.009 

1.11 ± 0.04 

0.74 ± 0.04 

St. dev. 

+0.022 ± 0.008 

0.0828 ± 0.0441 

0.13 + 0.04 


where po is a density normalisation, A is the factor by which 
the halo is overdense with respect to the critical density of 
the universe at the object redshift, r s = 1 /ca is a charac¬ 
teristic scale length measured in units of ta, 7 is the inner 
slope, /3 is the outer slope and a = 1. All parameters are 
calculated in terms of A = 500. 


There is a strong degeneracy between /? and r 3 when 
fitting the generalised NFW profile, so we estimated a value 
of C500 for each cluster using equation ( 12 ) from Dolag et al. 
( 2004 ) to be used in subsequent calculations. The values 
are shown in Table [9] When analysing stacked profiles, we 
assume the mean value of C500 from all the contributing 
clusters. 
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Table 9. Concentration parameters ca with respect to the critical density. The values of C 200 ,m were calculated using equation (12) in 
|Dolag et al.| ( |2004| . C 200 and C 500 were derived from C200,m assuming our fiducial cosmology. 


Object 

C200,m 

C200 


C 500 

RXC J0006.0-3443 

7.46 ±0.07 

4.43 ± 

0.04 

2.907 ± 0.030 

RXC J0049.4-2931 

8.16 ±0.06 

4.87 ± 

0.04 

3.217 ± 0.028 

RXC J0345.7-4112 

9.04 ± 0.07 

5.43 ± 

0.04 

3.609 ±0.029 

RXC J0547.6-3152 

7.04 ± 0.07 

4.17 ± 

0.04 

2.724 ± 0.031 

RXC J0605.8-3518 

7.29 ± 0.07 

4.32 ± 

0.04 

2.833 ± 0.030 

RXC J0616.8-4748 

7.72 ±0.06 

4.60 ± 

0.04 

3.024 ± 0.029 

RXC J0645.4-5413 

6.67 ± 0.07 

3.93 ± 

0.05 

2.560 ± 0.032 

RXC J0821.8+0112 

8.57 ± 0.06 

5.14 ± 

0.04 

3.400 ± 0.028 

RXC J2023.0-2056 

8.89 ± 0.07 

5.33 ± 

0.04 

3.539 ±0.029 

RXC J2048.1-1750 

7.16 ±0.07 

4.24 ± 

0.04 

2.775 ± 0.030 

RXC J2129.8-5048 

8.15 ±0.07 

4.87 ± 

0.04 

3.213 ±0.029 

RXC J2218.6-3853 

7.10 ±0.07 

4.21 ± 

0.04 

2.751 ± 0.031 

RXC J2234.5-3744 

6.76 ±0.07 

3.99 ± 

0.05 

2.600 ± 0.033 

RXC J2319.6-7313 

8.31 ± 0.06 

4.97 ± 

0.04 

3.284 ± 0.028 

Mean 

7.737 ± 0.018 

4.608 ±0.011 

3.031 ± 0.008 

St. dev. 

0.785 ±0.019 

0.496 ±0.012 

0.346 ± 0.008 


3.7 Brightest cluster galaxy properties and 
positions 

Large cD galaxies in the cluster central region are expected 
to have a significant impact on the distribution of fainter 
galaxies, which are less tightly bound and more susceptible 
to disruption than larger bright galaxies. Their distances 
from the cluster centres and sizes are given in Table m 

The BCGs in REXCESS were studied by |Haarsma et ah] 


1 2010). In this paper we use variable aperture elliptical mag¬ 


nitudes, whereas |Haarsma et al.| used magnitudes measured 
in a 12 /i _1 kpc fixed aperture metric radius, and as a con¬ 
sequence the magnitudes we measure here are higher. 


3.8 Radial density profiles 


By measuring the shape of the radial density profiles of 
galaxy clusters, we can assess the extent to which galax¬ 
ies (effectively collisionless particles) trace gas (a collisional 
fluid) and dark matter (which we assume to be a collisionless 
fluid) in galaxy clusters of different morphological types and 
masses. Since the evolutionary history of a galaxy may be 
quite different if it is on the red-sequence as opposed to in 
the blue cloud, or bright as opposed to faint, we must con¬ 
sider each sub-population of galaxies (red and bright, red 
and faint, blue and bright, and blue and faint) separately. 

We produced radial surface density profiles S for in¬ 
dividual clusters by summing object counts in radial bins 


around the cluster centres (see Section {'2.2 for details on 
the centring method), with annular bin edges at ( 0 , 0 . 1 , 
0.15, 0.2, 0.25, 0.35, 0.5, 0.7, 1, 1.25, 1.5, 2.5 and 5) r 5 oo- 
The number counts were divided by the total exposed area 
in each annulus to get the surface density. 


3.8.1 Count density profile fitting 

We fitted the generalised NFW model profile defined in 
Equation |5| to each of the individual galaxy count density 
profiles. This measurement is used to constrain the back¬ 
ground count density and to provide an estimate of the total 
galaxy count within a particular radius. 


When fitting the model profile to projected quantities 
like the observed galaxy count density profiles, we integ¬ 
rate p(r) numerically along the line of sight, and allow an 
additional uniform surface density component 5b g NFW as 
background, which is not projected. We fixed a = 7 = 1 be¬ 
cause the statistics of the individual profiles were too poor 
to constrain these parameters. 

An additional limit, p(r > r cu t 0 ff) = 0 is imposed so 
that the integrals converge and to break the degeneracy 
between ft and S'bgNFW- We found that values of SbgNFW 
were consistent within their uncertainties for r cu toff = (2.5, 
5 and 10) rsoo for the red populations which have good stat¬ 
istics, but that the poor statistics of the blue populations 
led to some cases of negative S'bgNFW if r C utoff > 2.5 rsoo• 
Whilst there are indications that the cluster extends bey¬ 
ond 2.5 rsoo, our data outside that radius are too sparse to 
make reasonable measurements, so we set r cu t 0 ff = 2.5 rsoo■ 
This limit is also used where we estimate the projected total 
mass profile within the cluster, and the projected gas density 
profile. 

To generate estimates of the total number of galaxies n g 
in each cluster, we integrated the fitted density profiles out 
to the relevant va- This was done for each galaxy population 
(combinations of red/blue and bright/faint) as well as for the 
combined total population (red and blue, down to the faint 
limit). 

Count density profiles for all of the individual clusters 
are shown in Figures B1 to B14 (online material) and their 


best fitting radial density profile parameters are given in 
Table [Tl] (full length table is given in Table |B 1 | online ma¬ 
terial) . 


3.8.2 Background count density analysis 

Before stacking the count density profiles of the clusters, the 
density of galaxies which are in the same line of sight as the 
cluster but not physically bound to it must be estimated 
and subtracted. Having already established that the objects 
may be embedded in more dense regions of the large scale 
structure (in Section { 3.3.1 1 , we made a number of differ¬ 
ent estimates of the local background for each observation 
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Table 10. BCG parameters. i?ABabsK is the absolute magnitude of the BCG and Lr is the corresponding luminosity, r is the distance 
of the BCG from the X-ray peak. A c is the ratio of the area of the BCG to 7r(0.1 r^oo) 2 - «K r0 n is the semi-major axis of the elliptical 
aperture used to measure the BCG luminosity. 


Object 

RAB abs K 
mag 

l r 

10 11 L 0 

r 

^500 

A c 

^Kron 

kpc 

RXC J0006.0-3443 

-23.393 ± 0.025 

1.65 + 0.04 

0.0033 

0.11267+ 0.00008 

55.590 + 0.028 

RXC J0049.4-2931 

-23.019 + 0.026 

1.168 + 0.028 

0.0056 

0.06321 ± 0.00008 

21.180 + 0.020 

RXC J0345.7-4112 

-22.92 + 0.05 

1.07 + 0.05 

0.0017 

0.099439 ± 0.000031 

23.026 + 0.005 

RXC J0547.6-3152 

-23.124 + 0.023 

1.287 + 0.027 

0.0241 

0.05545 ± 0.00009 

28.507 + 0.034 

RXC J0605.8-3518 

-23.455 + 0.022 

1.745 ± 0.035 

0.0049 

0.08959 ± 0.00008 

35.869 + 0.023 

RXC J0616.8-4748 

-23.791 ± 0.023 

2.38 + 0.05 

0.0068 

0.14889+ 0.00010 

57.391 ± 0.028 

RXC J0645.4-5413 

-23.576 ± 0.024 

1.95 + 0.04 

0.0095 

0.10326+ 0.00011 

53.77 + 0.04 

RXC J0821.8+0112 

-22.625 + 0.025 

0.813 + 0.019 

0.0057 

0.042928 ± 0.000019 

18.431 ± 0.005 

RXC J2023.0-2056 

-22.204 + 0.028 

0.551 ± 0.014 

0.0093 

0.064183 ± 0.000027 

19.517 + 0.006 

RXC J2048.1-1750 

-23.777 + 0.023 

2.35 + 0.05 

0.1729 

0.07459 ± 0.00005 

30.889 + 0.016 

RXC J2129.8-5048 

-22.601 ± 0.024 

0.795 ± 0.018 

0.1130 

0.041114 + 0.000023 

21.526 + 0.009 

RXC J2218.6-3853 

-23.424 + 0.023 

1.70 + 0.04 

0.0283 

0.06235 ± 0.00013 

32.02 ± 0.05 

RXC J2234.5-3744 

-23.218 + 0.023 

1.404 + 0.030 

0.1475 

0.041763 ± 0.000029 

29.832 + 0.014 

RXC J2319.6-7313 

-22.127 + 0.031 

0.514 + 0.015 

0.0061 

0.06545 ± 0.00027 

24.68 + 0.07 


Table 11. NFW fitting results for each of the galaxy population filters - bright and faint; red sequence and blue. The full length table 
is given in Table |B1| (online material). 


Object 

Galaxy filter 

p 

So 

r ~ 3 

'500 

SbgNFW 

r“ 2 
' 500 

RXC J0006.0-3443 

Bright, red 

2.7 + 0.5 

49.7 + 9.0 

7.2 + 2.1 

RXC J0049.4-2931 

Bright, red 

3.4 + 0.5 

59.0+ 10.9 

1.9 ±0.5 

RXC J0345.7-4112 

Bright, red 

3.2+ 1.0 

12.2 + 4.5 

2.3 ±0.6 

RXC J0547.6-3152 

Bright, red 

2.67 + 0.35 

50.1 ± 6.8 

5.4+ 1.4 

RXC J0605.8-3518 

Bright, red 

3.2 ±0.4 

55.5 + 8.8 

7.2 ±0.8 

RXC J0616.8-4748 

Bright, red 

2.58 + 0.26 

30.1 ± 3.4 

2.7 + 0.5 

RXC J0645.4-5413 

Bright, red 

2.74 + 0.29 

96.2+ 13.0 

9.8+ 1.5 

RXC J0821.8+0112 

Bright, red 

2.5 ±0.4 

25.8 + 4.0 

1.8 ±0.9 


set, with the aim of assessing whether the background level 
was high because the cluster itself is extended, or due to 
a uniform surface density of objects across the observation 
field. 

The simplest estimate Sbgsimpie was found by taking the 
count density in the region bounded by r > 1.5rsoo. This 
method is susceptible to contamination by the wings of the 
galaxy density distribution, including virialised galaxies and 
structures like infalling subclusters and filaments. 

We control for infalling structures which are not azi- 
muthally symmetric by making n independent estimates 
•Sbgsector (with standard deviation crs b ctor ) of the back¬ 
ground level in annular sectors in the same bounded re¬ 
gion. If the objects in the sectors are spatially uncorrel¬ 
ated and obey Poisson statistics, we would expect that the 
combined uncertainty of the n independent samples of the 
background density, a Sbgaectot = cT Sbgaecto Jx/n, would be 
the same as the Poissonian value for the complete region, 
i.e. as h g8 ec t0 r ~ Q Sb g si mp i 0 - O n the other hand, spatially cor¬ 
related objects are not drawn from a single Poisson distri¬ 
bution and appear significantly more frequently in some the 
sectors, leading to a higher Qs b soctor , i.e. we would expect 
O' S'bg sect or > os bgsimple if local structures are present. 

To control for broad wings of the galaxy density dis¬ 
tribution, we fitted the NFW model described in Equa¬ 
tion |5} - including a constant surface density component 
SbgNFW “ to the individual cluster profiles. SbgNFW is sens¬ 
itive to uniform density across the field, whereas Sbgsimpie 
includes the cluster wings and infalling structures as well. 


Sbgsimpie/‘S'bgNFW + 1 indicates that ‘Sbgsimpie is contam¬ 
inated by the wings of the cluster and overestimates the 
background level. 

The results of these measurements for the red galaxies 
down to the faint magnitude limit are given in Table [T2{ and 
Table [Pf] includes comprehensive results for all of the galaxy 
populations used in the radial profiles (the full length table 
is given in Table Cl online material). 

Comparing Qs bg sector and Q s bgS impie! we see that 7 of 
the clusters have aSbgssctor/ a Sbgsim P ie > 1-5, consistent with 
the existence of localised substructures within some of the 
sectors. We conclude that in half of the sample, filamentary 
structure or infalling objects contribute to the high back¬ 
ground in the cluster outskirts. 

Values of > 1 for the ratio Sbgahnpie/SbgNFW are sug¬ 
gestive that, in most cases, Sbgsimpie is an overestimate con¬ 
taminated by the wings of the cluster profile. (In the case of 
RXC J2023.0-2056, SbgNFW is consistent with zero.) Meas¬ 
urements of the remnant luminosity function in this region, 
once a model of the field galaxy density is removed, are 
presented in Section §3.9.2| 

We can also compare Sb g sim P ie/SbgNFW with f s 


from Table [5] In RXC J0345.7-4112, RXC J0645.4- 
5413, RXC J0821.8+0112 and RXC J2129.8-5048, 
Sbgsimpie/SbgNF^v f &. so at least a portion of the 

overdensity seen in the number counts function is likely due 
to contamination from the cluster itself, rather than the 
large scale structure along the line of sight. 

The least contaminated estimate of the background 
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Table 12. Background count density measurements for the red bright and faint galaxy population: simple count density Sbgsimple 
including Poisson uncertainties; standard error measured from sectors as bgS e ctori constant density component from an NFW fit SbgNFW- 
Comprehensive results for red-bright, red-faint, blue-bright and blue-faint galaxy populations are given in Table [T3| 


Object 

Sbgsimple 

arcmin -2 

^Sbgsector 

arcmin -2 

SbgNFW 

arcmin -2 

Sbgsimple/SbgNFW 

^Sbgsector /^Sbgsimple 

RXC J0006.0-3443 

0.86 ±0.05 

0.04 

0.74 ± 0.08 

1.17 

0.9 

RXC J0049.4-2931 

0.490 ± 0.029 

0.04 

0.471 ± 0.019 

1.04 

1.4 

RXC J0345.7-4112 

0.124 ± 0.021 

0.02 

0.090 ± 0.035 

1.38 

1.1 

RXC J0547.6-3152 

0.89 ± 0.05 

0.08 

0.78 ±0.08 

1.15 

1.6 

RXC J0605.8-3518 

1.11 ± 0.04 

0.08 

1.07 ± 0.06 

1.04 

1.9 

RXC J0616.8-4748 

0.73 ±0.04 

0.04 

0.67 ± 0.06 

1.09 

1.1 

RXC J0645.4-5413 

1.49 ± 0.05 

0.13 

1.30 ± 0.10 

1.15 

2.7 

RXC J0821.8+0112 

0.356 ±0.029 

0.05 

0.28 ±0.05 

1.27 

1.7 

RXC J2023.0-2056 

0.13 ±0.04 

0.04 

0.05 ±0.10 

2.41 

1.0 

RXC J2048.1-1750 

1.61 ± 0.06 

0.10 

1.38 ± 0.16 

1.17 

1.8 

RXC J2129.8-5048 

0.34 ± 0.04 

0.07 

0.284 ± 0.031 

1.20 

1.8 

RXC J2218.6-3853 

0.80 ± 0.04 

0.07 

0.81 ± 0.05 

0.99 

1.8 

RXC J2234.5-3744 

1.15 ±0.05 

0.05 

1.02 ±0.06 

1.13 

1.1 

RXC J2319.6-7313 

0.519 ±0.030 

0.02 

0.47 ± 0.05 

1.11 

0.7 


Table 13. Background count density measurements. In the case of the RXC J2023.0-2056 bright blue filter, no objects are detected in 
the region used for measuring Sbgsimple- The full length table is provided in Table |Cl| (online material). 


Object 

Galaxy filter 

Sbgsimple 
arcmin -2 

^Sbgsector 

arcmin -2 

SbgNFW 
arcmin -2 

Sbgsimple/SbgNFW 

RXC J0006.0-3443 

Bright, red 

0.166 ±0.020 

0.02 

0.12 ± 0.04 

1.34 

RXC J0006.0-3443 

Faint, red 

0.70 ± 0.04 

0.04 

0.62 ± 0.06 

1.12 

RXC J0006.0-3443 

Bright, blue 

0.128 ±0.017 

0.02 

0.102 ±0.029 

1.25 

RXC J0006.0-3443 

Faint, blue 

1.14 ±0.05 

0.07 

0.66 ±0.18 

1.73 

RXC J0049.4-2931 

Bright, red 

0.064 ± 0.011 

0.01 

0.051 ± 0.013 

1.25 

RXC J0049.4-2931 

Faint, red 

0.427 ± 0.027 

0.03 

0.411 ± 0.019 

1.04 

RXC J0049.4-2931 

Bright, blue 

0.074 ± 0.011 

0.01 

0.054 ± 0.022 

1.38 

RXC J0049.4-2931 

Faint, blue 

0.81 ± 0.04 

0.04 

0.795 ±0.024 

1.02 


galaxy count density appears to be SbgNFW i so we use that 
to generate background subtracted count density profiles. 
Whilst relying on a background measurement which is de¬ 
pendent on a prior cluster shape rather than one which is 
shape independent is not ideal, there is sufficient scope in 
the NFW model to make a reasonable characterisation of 
most of the clusters we see. 


the same r cuto ff. Typically, the projected electron density 
profiles are dominated outside of ~ 0.8 rsoo by extrapolated 
densities. 

Since the electron, total mass and galaxy count dens¬ 
ity profiles each have different normalisation, we plot 
S/Y, ^ S a, where a is the total area in the annulus, 
wherever they are compared. 


3.8.3 Comparison profiles and normalisation 

With the given count statistics and binning, a deprojection 
of the count density profiles produces very noisy profiles and 
therefore comparisons are best made with other projected 
profiles. 

The deprojected electron density profiles of |Croston 


et al. (2008) do not extend to r cuto fj (the edge of our galaxy 
density model) because the cluster X-ray emission becomes 
too faint to observe, so these were extrapolated using a gen¬ 
eralised NFW profile (described in Equation <§) fit to Cro- 
ston’s data, assuming C500 from Table [9] and allowing /3, 7 
and po to vary. 

We also wished to compare the galaxy count density 
profiles with some model of the total mass of the cluster. 
Since this is dominated by dark matter, we assumed the 
NFW density profile form described in Equation with 
ft = 3 and a = 7 = 1 , and the relevant C500 from Table [ 9 ] 

The electron and total mass profiles were projected us¬ 
ing the same numerical integrator as used for the galaxy 
count density profile fitting, in the same radial bins and with 


3.8.4 Stacked profiles 

We produce stacked galaxy density profiles for all clusters 
and each subsample by taking the mean of Snorm = 
(S — SbgNFw) / n g,r+h where n g , r +b is the total galaxy count 
for the cluster measured from the radial profile of all red 


and blue galaxies down the faint limit (see Section 6 3.8.1). 
The units of S'norm and S are r^ 0 . We produce stacked 
electron density profiles by taking the mean of S e ,norm = 
S e / J2r<r 500 S <? where S e ,norm and Se have units r^ 0 2 0 . Sim¬ 
ilarly, we produce stacked total mass profiles by taking the 
mean of Sm, norm — Sm / ^ „ f| Sm a , where Sat, norm has 

units rg 00 and Sm has units Mq r$ 00 . The mean profiles for 
the whole sample are shown in Figure [ 6 ] 

The bright red profile follows the total mass profile reas¬ 
onably well. The faint red profile is a little broader and there 
is a 3a (3x the galaxy count density uncertainty) deficit of 
faint red galaxies in the central bin. The blue profiles are 
both substantially broader than the total mass profile. There 
appears to be suppression of faint blue objects compared to 
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Figure 6. Mean projected radial profile of all of the clusters. Galaxy number density profiles are shown as black points with error bars, 
with the best fitting NFW model shown as a red dotted line. The projected electron density is shown as a blue dashed line. The NFW 
profile with parameters (a = 1,/3 = 3 ,7 = 1) and C 500 from Table [ 9 ] is shown in grey. The upper panels represent galaxies on the red 
sequence, and the lower panels include only galaxies bluer than the red sequence. The panels on the left are ‘bright’ and the panels on 
the right are ‘faint’ galaxies. The radial densities S shown are unitless, having been normalised by the mean S of all of the annuli within 
7*500, such that they overlay each other as closely as possible. 


the overall trend which is limited to the region < 0.15 7 * 500 ; 
with respect to the NFW total mass curve (shown in grey) 
the deficit is around 5cr. The presence of a cD brightest 
cluster galaxy in the cluster core could make the detection 
of faint objects more difficult, but since the BCG occupies a 
small fraction of the area of the central bin - if indeed it is 
positioned there - this effect is negligible (see Table [To| for 
BCG sizes relative to the central bin area, and the distance 
of the BCG from the X-ray peak). Additionally, we would 
expect to see a similar effect for the red sequence galaxies, 
for which there is no evidence. 

I 11 Section §2.4| we described subsamples of the REX¬ 
CESS dataset based on morphological classifications and 
and cluster total mass. Mean profiles for the massive and low 
mass subsamples are shown in Figure [7] for the disturbed 
and regular subsamples in Figure [8] and for the cool-core 
and non-cool-core subsamples in Figure [9] 

The subsample profiles appear to show that the sup¬ 
pression of faint blue galaxies in the cluster cores is driven 
by the massive subsample where this effect is marked, and 
is absent in the low mass subsample (both in Figure[7|. The 
massive subsample also shows suppression of the faint red 
galaxies in the innermost radial bin (< O.lrsoo), not seen 
in the mean profile for the whole sample. Figure [9] shows 
that there is even stronger suppression of faint blue galax¬ 
ies in the cores of the cool-core clusters (only one of which 


is classified as massive), and this effect appears in all three 
cool-core clusters’ individual profiles (see Section |B] online 
material). Additionally, suppression of bright blue galaxies 
is noted in the cool-core cluster cores, an effect not seen at 
all in the stacked profile for the full sample. In the cool-core 
clusters the red populations appear to be unaffected. 

There is some suppression of faint blue galaxies in the 
centres of the regular clusters, but other than that the pro¬ 
files for the regular and disturbed clusters appear qualit¬ 
atively very similar to one another. In particular, we find 
no evidence from this analysis that there is any substantial 
difference between the two subsamples which might give a 
measure of the dynamical state to complement the X-ray 
based centre shifts parameter used for the disturbed/regular 
classification (see Section j |2.4[ ). The difference which is seen 
could be statistical noise. 

There is strong evidence that the profiles extend at least 
up to the 2.5 r 500 limit, as the outer bins of the stacked 
profiles have a significant positive residual even after back¬ 
ground subtraction. 

We fitted projected generalised NFW models to the 
stacked galaxy count density profiles using the same method 
as used for the individual clusters to yield /3 g , allowing for a 
constant density component SbgNFW in case the background 
subtraction before stacking was incomplete. We include for 
comparison the mean electron density profile outer slope /3 e , 
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(a) Massive clusters: RXC J0006.0-3443, RXC J0547.6-3152, 
RXC J0605.8-3518, RXC J0645.4-5413, RXC J2048.1-1750, 
RXC J2218.6-3853, RXC J2234.5-3744 


(b) Low mass clusters: RXC J0049.4—2931, RXCJ0345.7— 
4112, RXC J0616.8-4748, RXC J0821.8+0112, RXC J2023.0- 
2056, RXC J2129.8-5048, RXC J2319.6-7313 


Figure 7. Stacked, projected radial number density profiles for clusters above and below the median mass in the population. The lines 
are the same as described in Figure [6] 




(a) Disturbed clusters: RXC J0006.0—3443, RXCJ0616.8- 
4748, RXC J2023.0-2056, RXC J2048.1-1750, RXCJ2129.8- 
5048, RXC J2218.6-3853, RXC J2319.6-7313 


(b) Regular clusters: RXC J0049.4-2931. RXC J0345.7—4112, 
RXC J0547.6-3152, RXC J0605.8-3518, RXC J0645.4-5413, 
RXC J0821.8+0112, RXC J2234.5-3744 


Figure 8. Stacked, projected radial number density profiles for disturbed and regular clusters. The lines are the same as described in 
Figure |6| 
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(a) Cool-core clusters: RXC J0345.7—4112, RXC J0605.8-3518, 
RXC J2319.6-7313 


(b) Non-cool-core clusters: RXC J0006.0-3443, 
2931, RXC J0547.6-3152, RXC J0616.8-4748, 
5413, RXC J0821.8+0112, RXC J2023.0-2056, 
1750, RXC J2129.8-5048, RXC J2218.6-3853, 
3744 


RXC J0049.4- 
RXC J0645.4- 
RXC J2048.1- 
RXC J2234.5- 


Figure 9. Stacked, projected radial number density profiles for cool-core and non-cool-core clusters. The lines are the same as described 
in Figure [6] 


which is calculated by taking the mean f3 e of all the clusters 
in the subsample. The parameters are tabulated in Table p~4] 
and the fits are shown in Figures [5] to [9] 


The best fitting profiles for both bright and faint red 
sequence galaxies have outer slopes which are flatter than, 
but nevertheless in rough agreement with j3 = 3, the slope 
of the assumed total mass profile. The bright blue profile is 
substantially broader than the NFW, and the best fitting has 
outer slope inconsistent with (5 = 3. The faint blue profile 
is similar, with a best fitting outer slope consistent with 
the outer slope for the bright blue galaxies but with large 
uncertainties. 


The cumulative fraction of red galaxies for the full 
sample, and for the high and low mass clusters is shown 
in Figure |10| These measurements reflect the morphology- 


density relation for ellipticals and spirals (e.g. Dressier 


19841, and our measured blue fraction at the limit of our 


observations - well within the cluster region of influence - is 
substantially lower than the field spiral population. Compar¬ 
ing the low mass and massive clusters, we see that outside 
0.2 r 500 , the red fraction reaches a plateau in the low mass 
clusters, but in the massive clusters it is still higher than 
90% and doesn’t reach the same plateau level even at the 
limit of our observations. Even in low mass clusters, the red 
fraction approaches 100% in the central regions. 


3.9 Luminosity measurements 

3.9.1 Variation of the cluster luminosity function with 
radius 

Since we see a reduction in the faint galaxy counts in 
cluster centres, we produced background subtracted lumin¬ 
osity functions for the projected annuli with edges at (0, 
0.15, 0.5 and 1.0) rsoo for the full galaxy population and for 
the red sequence galaxies. A selection of these luminosity 
functions are shown in Figure[lT] The method of generation, 
normalisation and stacking is described in Section §3.3.2| the 
only difference is that we now impose an additional cata¬ 
logue selection based on the red sequence fit. 

Figure [lT] shows that the red sequence luminosity func¬ 
tions outside of 0.15 rsoo are all extremely similar. The in¬ 
ner luminosity functions have a break at around Rab abs k = 
— 18, above which the function drops below the trend. This 
suppression is largely due to the massive clusters. The un¬ 
certainties on the inner luminosity function for the low mass 
clusters are too large to conclude that there is suppression; 
within the uncertainties it appears that the faint galaxy 
count continues the trend seen at brighter magnitudes. We 
note that the massive clusters are more distant on average, 
and that the magnitude limit is lower for these observa¬ 
tions, but not sufficiently low that it explains the break at 
Rab abs k = —18. Additionally, the stacking procedure - nor¬ 
malising to a complete part of the luminosity function and 
ensuring that truncated magnitude bins do not contribute to 
the mean - should minimize any influence of the complete¬ 
ness limit on the final luminosity function shapes. There is 
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Rab abs K /mag 


(a) Stacked luminosity functions for all the clusters in the data- (b) Stacked luminosity functions for for all the clusters in the 

set, including red sequence galaxies satisfying |iu ra | < 3. dataset, including blue galaxies satisfying ut ra < —3. 




(c) Stacked luminosity functions for massive clusters, including (d) Stacked luminosity functions for low mass clusters, includ- 

red sequence galaxies satisfying |u) r s| < 3. ing red sequence galaxies satisfying \w T s < 3. 

Figure 11. Stacked luminosity functions for the annuli bounded by (0, 0.15, 0.5 and 1.0)r5oo- The functions are normalised and then 
artificially separated by a factor of 10. 


some evidence that the massive sample luminosity function 
is a little flatter than the low mass sample one, but given 
the size of the uncertainties it is difficult to be certain. 

The blue object luminosity functions vary strongly with 
respect to cluster-centric distance. Figure [l lb| shows an ex¬ 
cess of bright galaxies in the outer cluster regions (0.5 < 
r/rsoo < 1) which isn’t seen at smaller radii, and we found 
that there was an excess of bright blue galaxies in the off- 
target region as well. In the two inner regions sampled 
(f/p 500 < 0.15 and 0.15 < r/rsoo < 0.5) there is strong 


variation in shape of the luminosity function, away from a 
simple Schechter function. 


There is no evidence of a difference between the lumin¬ 
osity function of disturbed and regular clusters, shown in 
Figure [12] Any differences which are apparent are consist¬ 
ent with being statistical effects. 
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(a) Stacked luminosity functions for disturbed clusters, includ- (b) Stacked luminosity functions for regular clusters, including 

ing red sequence galaxies satisfying \w rs \ < 3. red sequence galaxies satisfying \w rs \ < 3. 

Figure 12. Stacked luminosity functions for the annuli bounded by (0, 0.15, 0.5 and 1.0) rsoo- The functions are normalised and then 
artificially separated by a factor of 10. 


Table 14. Stacked cluster best fitting parameters. 1 No error is 
quoted where the fitting routine failed to estimate the covariance 
matrix. 


Object Galaxy filter (3g (3 e 


All clusters 

Massive 

Low mass 

Disturbed 

Regular 

Cool-core 

Non-cool-core 

All clusters 

Massive 

Low mass 

Disturbed 

Regular 

Cool-core 

Non-cool-core 

All clusters 

Massive 

Low mass 

Disturbed 

Regular 

Cool-core 

Non-cool-core 

All clusters 

Massive 

Low mass 

Disturbed 

Regular 

Cool-core 

Non-cool-core 


Bright, red 
Bright, red 
Bright, red 
Bright, red 
Bright, red 
Bright, red 
Bright, red 
Faint, red 
Faint, red 
Faint, red 
Faint, red 
Faint, red 
Faint, red 
Faint, red 
Bright, blue 
Bright, blue 
Bright, blue 
Bright, blue 
Bright, blue 
Bright, blue 
Bright, blue 
Faint, blue 
Faint, blue 
Faint, blue 
Faint, blue 
Faint, blue 
Faint, blue 
Faint, blue 


2.78 ±0.16 
2.97 ±0.16 
2.57 ±0.30 
2.53 ±0.24 
3.05 ±0.13 

2.3 ±0.4 
2.86 ±0.17 
2.51 ±0.23 

2.4 ±0.4 
2.55 ±0.32 
2.29 ±0.28 
2.68 ±0.24 
2.6 ±0.5 
2.40 ±0.25 

1.79 ±0.26 
- 0.1 ± 1.6 
1.84 ±0.20 
1.63 ±0.34 
2.07 ±0.34 
1.29 1 

1.4 ±0.5 

1.3 ±0.4 

1.4 ± 0.6 
0.6 ± 1.2 
0.7 ± 0.6 
1.8 ±0.6 
1.2 ± 1.8 
1.14 ± 0.35 


2.936 ±0.026 
3.131 ± 0.032 
2.74 ± 0.04 
2.853 ± 0.034 
3.02 ± 0.04 
2.70 ± 0.06 
3.001 ± 0.028 
2.936 ±0.026 
3.131 ± 0.032 
2.74 ± 0.04 
2.853 ±0.034 
3.02 ±0.04 
2.70 ±0.06 
3.001 ± 0.028 
2.936 ±0.026 
3.131 ± 0.032 
2.74 ± 0.04 
2.853 ±0.034 
3.02 ± 0.04 
2.70 ± 0.06 
3.001 ± 0.028 
2.936 ±0.026 
3.131 ± 0.032 
2.74 ± 0.04 
2.853 ± 0.034 
3.02 ± 0.04 
2.70 ± 0.06 
3.001 ± 0.028 


3.9.2 Residual cluster luminosity function in the off-target 
region 

Given the evidence of structures in the off-target region in 
Section §3.8.2| and Section §3.8.4[ we re-analysed the lumin¬ 
osity function in the r 2 oo < r region to try to find a residual 
cluster luminosity function, once our assumed field galaxy 
function was subtracted. 

The possible excess of cluster galaxies outside rsoo 
means that the normalisation factor / g found for the as¬ 


sumed Held galaxy function ^Metcalfe (see Section £ 3.3.11 may 
be slightly overestimated. However, since the cluster lumin¬ 
osity function is largely invariant with distance from the 
cluster centre, we can attempt to fit both cluster and field 
simultaneously for the full image. This should give an im¬ 
proved estimate of / g , which can be used with ^Metcalfe as 
the background. 

Assuming values for a and M , from Table [ij] we fit¬ 
ted a combined model £ fa n off (<f>(L) ± / g ^Metcalfe) to the lu¬ 
minosity histogram of each full field, making no magnitude 
or red-sequence based selections since these may alter the 
field number counts. The ^Metcalfe component was subtrac¬ 
ted from the count histogram in the off-target region, and 
the results are shown in Figure ^3] 

The shape of the residual in Figure [l3]appears, in most 
cases, inconsistent with the shape of the background num- 


4 a and M* were fitted to the entire sample. Since the cluster 
luminosity function is invariant with distance from the cluster 
centre, the background subtraction performed before the stacking 
operation in Section §3.3.2| should not bias the shape of the final 
luminosity function, despite the small residual of cluster galaxies 
in the region used as the background. 
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Figure 13. Off-target region (r200 < r ) residual cluster luminosity function. The assumed background model (including the falloff 
component) is shown as a dashed line for comparison. 


ber counts function which is also shown. The residual often 
resembles the cluster luminosity function. This is consist¬ 
ent with the tentative conclusion drawn in Section §3.8.2[ 
that the cluster does extend some distance outside of r 2 oo- 
Compared with the on-target luminosity functions shown in 
Section §3.9.1[ these residuals have denser faint components, 
with an upturn at fainter magnitudes similar to the dwarf 
upturn seen in Popesso et al. (2006), but are also consistent 
with being due to some remnant field contamination. 


3.9.3 Total cluster luminosity 

The total cluster luminosity and the related mass-to-light 
ratio are useful parameters when assessing the efficiency 
or disruption of star formation in different types of galaxy 
clusters. Since we are dealing with projected data, a cor¬ 
rection needed to be made for the cluster galaxies outside 
of rsoo or r 2 oo but, when seen in projection, in one of the 
annular bins. 

We considered two populations of galaxies when cal¬ 
culating the total cluster luminosities: |w ra | < 3 which in¬ 
cludes only galaxies on the red sequence, and w IS < 3 which 
includes both the red sequence and the blue cloud, but ex¬ 
cludes objects redder than the red sequence. To ensure that 
the correction to the total luminosity required due to the 
magnitude limit was approximately equal for all the clusters 
in the sample, we considered galaxies satisfying R < R, + 5. 

Using the best fitting NFW models - which include a 
background estimate - for each of the galaxy populations, 
we assign a weight w which represents the probability that 


a galaxy seen in a particular radial bin is within ta- If the 
total best fitting model count density for a particular radial 

bin is fft otai = ■S'bgNFW + iS'NFW,r cutof f=2.5 t - 500 (as described in 

Section 53.8.11, the weight assigned to each galaxy in that 
radial bin is wa = 5NFW,r cutoff =r A /Stotai, where Snfw is 
found by integrating the volume density model p in annuli 
along the line of sight, and setting p = 0 where r > r cu toff- 
The total luminosity within va is then the sum JT wa i Li, 
and the total count is JT w i f° r a ll galaxies i in the par¬ 
ticular population. The major source of uncertainty in this 
calculation is the uncertainty on w, arising from the uncer¬ 
tainties on the fitting parameters in the model. The BCG is 
assigned w = 1, but other galaxies are not specially treated. 
Typical values of w are around 0.9 in the innermost radial 
bins. 


Since the radial count density profiles for the bright red, 
faint red, bright blue and faint blue populations are different, 
we calculate luminosities for all four subpopulations separ¬ 
ately, and then sum the relevant sub-populations to get total 
red or red plus blue luminosities. In cases where the best fit¬ 
ting model is consistent with there being no overdensity for a 
particular population and has very large uncertainties on the 
relevant fit parameters this subpopulation is not included in 
the final total luminosity. Only the RXC J0345.7-4112 faint 
blue population is affected by this step. 


The total luminosities for each cluster are given in 
Table H5l 
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Table 15. Total R band luminosities. 


Object 

Galaxy filter 




7-500 

10 12 L ( 

3 

7-200 

1O 12 L 0 

RXC J0006.0-3443 

Bright 

and 

faint, 

red 



1.49 ±0.16 

2.21 ± 0.25 

RXC J0006.0-3443 

Bright 

and 

faint, 

red 

and 

blue 

1.78 ±0.17 

2.83 ±0.28 

RXC J0049.4-2931 

Bright 

and 

faint, 

red 



1.02 ±0.12 

1.30 ±0.15 

RXC J0049.4-2931 

Bright 

and 

faint, 

red 

and 

blue 

1.23 ±0.16 

1.74 ±0.24 

RXC J0345.7-4112 

Bright 

and 

faint, 

red 



0.266 ±0.028 

0.33 ± 0.04 

RXC J0345.7-4112 

Bright 

and 

faint, 

red 

and 

blue 

0.31 ± 

0.04 

0.38 ±0.05 

RXC J0547.6-3152 

Bright 

and 

faint, 

red 



2.08 ± 

0.17 

2.87 ± 0.25 

RXC J0547.6-3152 

Bright 

and 

faint, 

red 

and 

blue 

2.20 ± 

0.22 

3.08 ±0.31 

RXC J0605.8-3518 

Bright 

and 

faint, 

red 



1.13 ± 

0.08 

1.46 ±0.11 

RXC J0605.8-3518 

Bright 

and 

faint, 

red 

and 

blue 

1.16 ± 

0.09 

1.51 ± 0.12 

RXC J0616.8-4748 

Bright 

and 

faint, 

red 



0.96 ± 

0.06 

1.43 ±0.10 

RXC J0616.8-4748 

Bright 

and 

faint, 

red 

and 

blue 

1.17 ± 

0.11 

1.79 ±0.17 

RXC J0645.4-5413 

Bright 

and 

faint, 

red 



3.59 ± 

0.26 

5.1 ± 0.4 

RXC J0645.4-5413 

Bright 

and 

faint, 

red 

and 

blue 

3.61 ± 

0.26 

5.1 ± 0.4 

RXC J0821.8+0112 

Bright 

and 

faint, 

red 



0.75 ± 

0.07 

1.18 ±0.12 

RXC J0821.8+0112 

Bright 

and 

faint, 

red 

and 

blue 

0.97 ± 

0.09 

1.63 ±0.17 

RXC J2023.0-2056 

Bright 

and 

faint, 

red 



0.49 ± 

0.11 

0.70 ±0.18 

RXC J2023.0-2056 

Bright 

and 

faint, 

red 

and 

blue 

0.63 ± 

0.15 

0.88 ±0.23 

RXC J2048.1-1750 

Bright 

and 

faint, 

red 



2.14 ± 

0.23 

3.4 ± 0.4 

RXC J2048.1-1750 

Bright 

and 

faint, 

red 

and 

blue 

2.63 ± 

0.24 

4.5 ±0.4 

RXC J2129.8-5048 

Bright 

and 

faint, 

red 



0.91 ± 

0.13 

1.14 ±0.19 

RXC J2129.8-5048 

Bright 

and 

faint, 

red 

and 

blue 

0.94 ± 

0.13 

1.22 ±0.19 

R.XCJ2218.6-3853 

Bright 

and 

faint, 

red 



1.04 ± 

0.21 

1.30 ±0.25 

R.XCJ2218.6-3853 

Bright 

and 

faint, 

red 

and 

blue 

1.19 ± 

0.21 

1.55 ±0.27 

RXC J2234.5-3744 

Bright 

and 

faint, 

red 



3.34 ± 

0.13 

4.38 ±0.17 

RXC J2234.5-3744 

Bright 

and 

faint, 

red 

and 

blue 

3.81 ± 

0.19 

5.4 ± 0.4 

R.XCJ2319.6-7313 

Bright 

and 

faint, 

red 



0.45 ± 

0.05 

0.77 ± 0.09 

RXC J2319.6-7313 

Bright 

and 

faint, 

red 

and 

blue 

0.56 ± 

0.07 

0.91 ± 0.11 


3.9.4 Mass-to-light ratio relation 

M 500 and M 200 were calculated using rsoo from Table [TJ con¬ 
centration parameters from Table [ 9 ] (to transform between 
rsoo and r 200 ) and the fiducial cosmology. We calculated 
mass-to-light ratios Ma/La for each of the clusters, which 

Using Ma/La = 


14 


are plotted against Ma in Figure 

7? (M A /M pivot r as a model, with M pivot = 5 x 1O 14 M 0 , we 
found best fitting parameters to the mass-to-light vs. mass 
relation, which are given in Table [I6| 

Three objects - RXC J0345.7-4112 (cool-core), 
RXC J00605.8-3518 (cool-core and massive) and 
RXC J2218.6-3853 (massive and disturbed) - lie slightly 
above the fitted mass-to-light relationship, but do not 
significantly affect the fit. 

We find a slope e of 0.06 ± 0.10 for the red sequence 
within 7*200, and 0.16 ± 0.14 for the red plus blue galaxy 
population within 7*200- The increase in slope when the blue 
luminosity is included compared with the case with just the 
red sequence luminosity is consistent with a decrease in blue 
fraction at high masses, already noted in Section §3.8.4| and 
shown in Figure [T0| The increase in the mean luminosity of 
blue galaxies as cluster mass increases and faint galaxies are 
disrupted/otherwise suppressed is insufficient to compensate 
for the decreased blue fraction. 

Popesso et al. ( |2007j ) measured M 200 /L 200 for red se¬ 
quence objects and quote a slope e = 0.18 ± 0.04, once pro¬ 
jection effects are taken into account. Whilst this result is 
in fair agreement with our measurements given the uncer¬ 
tainties, both of our best estimates for the red sequence are 
somewhat flatter. 


mass-to-light ratio of e = 0.33 ± 0.02 for objects in the 
MaxBCG catalogue of galaxy clusters in the Sloan Digital 
Sky Survey measured in the i band, a value in fair agreement 
with our measurement for the red and blue populations, but 
in this case too, our result is flatter. 

The M 200 range available to both [Popesso et al.| and 

A r. I .. 1 Jn -YX W — I A T 


Sheldon et al. 


Sheldon et al. (2009bI quote a logarithmic slope on the 


(5 x 10 12 - lO iJ )/i _1 M 0 , substantially 
larger than the single order of magnitude mass range in 
REXCESS, lending their analyses greater power to resolve 
mass dependent effects. 

Carlberg et al. (19971 quote an asymptotic value 
(289 ± 50) h M 0 Lq 1 for the Gunn r band, in excellent agree¬ 
ment with our value of (3.0 ± 0.4) x 10 2 h M 0 Lq 1 for the 
red and blue galaxies in the R band, measured at 10 15 M 0 . 


4 DISCUSSION 

Both the invariance of the cluster luminosity function with 
respect to radius (outside of cluster centres, i.e. r > 
0.15 7 * 500 ) and suppression of faint galaxies in the central 
regions of galaxy clusters have been noted before ( |Popesso| 
et al.p006 l. Popesso et al. found significant suppression in 
the late type luminosity function (corresponding to our blue 
population) for small cluster-centric distances; we find that 
there is substantial change in the red luminosity function 
close to the cluster core as well. 

The colour-magnitude relation parameters drawn from 
the WINGS clusters ( Valentinuzzi et al.|2011 I are quite sim¬ 
ilar to the values we see here, although the gradient scatter 
from our sample is twice as large as that from WINGS. That 
sample is also X-ray selected, but from clusters with lower 
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Figure 14. Mass-to-light relations. The marker labels are the IDs given in Table [T] The uncertainties only take into account the 
luminosities; mass uncertainties have not been considered. 


Table 16. Best fitting parameters to the mass-to-light ratio relation M A /L A = r\ (M A /M v \ vni y where M p i vot = 5 X 10 14 Mq. The 
values of the M A L~y 1 measured at 10 14 M@ and 10 15 Mq are corrected for the magnitude cut at R < 1?* + 5 using the correction factor 
described in Section §3.3.2| r/ is not corrected for the magnitude cut. 


Galaxy filter 

A 

V 

M 0 L-‘ 

e 

M a 1a (iO 14 ^) 
MoL- 1 

M a L^(10 15 Mq) 

M© L© 1 

M A L^ L h~ 1 (1O 15 M 0 ) 
MqL- 1 

Bright and faint, red 

200 

242.9 ± 16.7 

0.06 ±0.10 

(2.2 ±0.4) X 10 2 

(2.5 ±0.2) x 10 2 

(3.5 ±0.3) x 10 2 

Bright and faint, red 

500 

250.4 ±15.4 

-0.00 ± 0.09 

(2.4 ±0.4) x 10 2 

(2.4 ± 0.2) X 10 2 

(3.5 ±0.3) x 10 2 

Bright and faint, red and blue 

200 

192.5 ± 18.0 

0.16 ±0.14 

(1.4 ±0.4) x 10 2 

(2.1 ± 0.3) x 10 2 

(3.0 ±0.4) x 10 2 

Bright and faint, red and blue 

500 

216.2 ± 15.6 

0.10 ±0.10 

(1.8 ±0.3) x 10 2 

(2.3 ±0.2) x 10 2 

(3.2 ±0.3) x 10 2 


redshifts. Whilst it may be the case that there is a tight¬ 
ening of the distribution of red sequence parameters at red- 
shifts approaching 2 = 0, it is difficult to distinguish this 
effect from the increased measurement uncertainties intro¬ 
duced by increasing numbers of field galaxies in the same 
region of colour-magnitude space. 

Of particular interest to us was any indication that the 
galaxy density profile of disturbed clusters is also disturbed. 
Both the slope of the galaxy count density distribution (/9 g ) 
and the luminosity functions could have shown differences, 
but there is no significant evidence of a difference in either 
of these two properties in the disturbed and regular sub¬ 
samples. The similarity in the luminosity functions echoes 


the findings of De Propris et al. (20131 where luminosity 


functions of collisional and normal clusters in a sample se¬ 
lected by X-ray, optical and weak and strong lensing were 
studied. 

We suggest that two main processes can be invoked to 
explain the distribution of red and blue, bright and faint 
galaxies in clusters. Ram pressure stripping occurs as a 


galaxy moves with velocity v through the intracluster me¬ 
dium (ICM) with density p, and the gas in the galaxy is 


subjected to pressure P tx pv 2 (Gunn & Gott 1972b The 


pressure ablates cool gas from the halo, slowing star form¬ 
ation and turning blue galaxies redder. This effect should 
be more pronounced in regions of galaxy clusters with high 
gas densities, in particular in cool-cores. The galaxy infall 
velocity is related to the cluster mass M by v 2 oc M, so ram 
pressure stripping should also be stronger for more massive 
clusters. Because this process affects bright (as well as faint) 
galaxies, which dominate the total luminosity of the cluster, 
it should lead to a decrease in fraction of the cluster lumin¬ 
osity provided by the blue galaxy population as cluster mass 
increases. As it affects star formation as a whole, it should 
also lead to decreased overall star formation efficiency in 
more massive clusters and to a positive slope on the mass-to- 
light ratio relation measured using just red sequence galax¬ 
ies. 

The second main process, harassment, occurs as weakly 
bound galaxies interact tidally with more massive objects. 
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Figure 10. Cumulative fraction of red galaxies. The top two pan¬ 
els show (left) bright and faint galaxies for the whole sample and 
(right) only bright galaxies for the whole sample. The remaining 
panels show bright and faint galaxies for the specified subsamples. 


Parts of the weakly bound galaxy are stripped away, or the 
galaxy is completely disrupted. The remnants are a source of 
intracluster light (ICL). This process is strongest in regions 
where galaxy count densities and velocities are highest and 
affects more weakly bound (less massive/lower luminosity) 
galaxies more. Because the galaxy count densities in the 
central regions of the clusters are not strongly dependent on 
mass, this effect is expected to be less mass dependent than 
ram pressure stripping. 

There are several key pieces of evidence we can use 
to disentangle the two processes. The suppression of faint 
galaxies independent of the galaxy colour in the densest re¬ 
gions of the galaxy clusters, with the strongest effect in the 
most massive clusters, suggests harassment - a gravitational 
process independent of gas density and star formation in the 
affected galaxy - is important. The steeper mass-to-light ra¬ 
tio relation for the blue plus red galaxies vs. the red sequence 
alone, as well as decreasing blue galaxy fraction with higher 
mass, is evidence that ram pressure stripping - a process 
which primarily affects blue galaxies - is increasingly ef¬ 
fective in reducing star formation rates as infall velocities 
of galaxies rise. From the flatter blue galaxy count density 


profiles in all of the clusters, it is clear that the blue galaxy 
population does not survive long enough to relax into the 
cluster potential before being stripped of its cold gas and be¬ 
coming part of the red population. There is some evidence 
that the suppression of blue galaxies is most pronounced in 
the regions with the highest gas densities at the centres of 
cool-core clusters, but given the small sample size we cannot 
be certain that this is not a statistical anomaly. 

The REXCESS sample was selected by X-ray luminos¬ 
ity, ensuring that only clusters which are well evolved and 
have deep gravitational potential wells with hot, dense ICM 
are selected. This is in contrast to clusters in optically se¬ 
lected samples which are not always as well evolved, and 
consequently may not have a sufficiently dense ICM for ef¬ 
ficient ram pressure stripping. Bohringer et al. (20041 note 
that cool-cores in clusters are long-lived, which may allow 
more time for processes which disrupt galaxies and stop star 
formation from taking place. 

Both the red sequence and red plus blue mass-to-light 
ratio relation slopes we measure are flatter than in the liter¬ 


ature, compared to both X-ray selected samples (Popesso 


et al. 


2007I or optically selected samples (Sheldon et al. 


2009a). Given the scatter in the relation and the relatively 


large uncertainties on the best fitting parameters, as well as 
the fact that the REXCESS sample contains only clusters 
spanning one order of magnitude at the highest masses, it is 
impossible using these data to distinguish between the case 
where the differences between the slopes measured here and 
in the literature are purely statistical in nature, or due to 
different physical processes in the two samples - e.g. stronger 
ICM effects. 

The ICL has not been taken into account in this work, 
but if the relative density of the ICL in the centres of massive 
clusters were higher than in low mass clusters, then this 
would be further evidence for increased harassment. If we 
assume that 10% of the light of galaxy clusters is ICL (e.g. 


Zibetti et al. 2005), our mass-to-light ratio normalisations 
may be overestimated by a factor of ~ 1.1, leading to a cor¬ 
rection of comparable size to the normalisation uncertain¬ 
ties. However, based on the measurements of the BCG sizes 
and luminosities as compared with Haarsma et al. (20101 


described in Section §3.7[ it seems likely that a substantial 
fraction of the intracluster light is included in the BCG lu¬ 
minosity we measure, so the correction may well be smaller. 

[Gonzalez et al. (2007) and Gonzal ez et al.| (|2013) discuss 
the reduced efficiency of ICL generation in more massive 
objects which is coupled with a higher X-ray gas fraction. 
Zi betti et al. | ( [2005 ) find that the ICL surface brightness is 
correlated with BCG luminosity, but that the total fraction 
of light contributed by the ICL is almost independent of 
cluster richness and BCG luminosity. Given the open dis¬ 
cussion on the ICL light fraction as a function of cluster 
mass, it is too premature to include the effect of the ICL in 
the mass-to-light ratio in our results. 


5 SUMMARY AND CONCLUSIONS 

We have used a sample of 14 galaxy clusters from the REX¬ 
CESS survey to investigate radial density profiles of galaxies 
and intra-cluster medium. 

• The red galaxy density traces the dark matter dens- 
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ity closely outside of the cluster centres (in the region 
r > 0.15 rsoo)- The best fitting NFW model outer-slopes 
jSg are roughly consistent with /3 S = 3, with a best estimate 
/3 g = 2.78 ± 0.16, fitted to the stacked bright red sequence 
galaxy density profile of all the clusters. 

• The blue sequence count density profile slopes are sub¬ 
stantially shallower than the (3 = 3 total mass model, with 
a best estimate /3 g = 1.79±0.26, fitted to the stacked bright 
blue galaxy density profile of all the clusters. 

• The mean outer slope for the gas density profiles of the 
full sample is /3 e = 2.936 ± 0.026. Within the cluster centres 
the gas and dark matter profiles tend to diverge. 

• We find that faint blue galaxies are suppressed in the 
centres of massive and regular clusters. Faint red galaxies 
are also suppressed in the centres of massive clusters. Both 
bright and faint blue galaxies are heavily suppressed in the 
centres of cool-core clusters, but the faint red galaxies are 
unaffected. This is consistent with the idea that the suppres¬ 
sion of star formation is driven by ram pressure stripping of 
gas from galaxies, but that wholesale disruption of galaxies 
is caused by galaxy interactions in regions with high galaxy 
densities. 

• Our measurement of the logarithmic slope e of the 
galaxy cluster mass-to-light relation within r 2 oo of 0.16±0.14 

for all galaxies, measured in the R band, is in fair agree¬ 
ment with e = 0.33T0.02 from |Sheldon et aL] ( |2009~b| ) , meas¬ 
ured in the i band. Our measurement of the mass to light ra¬ 
tio normalisation of (3.0 ± 0.4) x 10 2 h Mg Lq 1 (evaluated 
at 10 15 Mg) in the R band is in excellent agreement with 
Carlberg et al. (1997|) measured in the Gunn r band. 


• There is no evidence of any difference in the galaxy 
count density profiles when comparing clusters classified as 
having disturbed X-ray morphology with those which are 
regular. 
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^Cluster/ r 500 ^Cluster/ r 500 


Figure Al. RXC J0006.0-3443 in the R band. Stars have been Figure A3. RXC J2234.5-3744 in the R band. Stars have been 
excised from this image. excised from this image. 


3 

6 



-10 1 
Cluster / r 500 



Cluster / r 500 


Figure A2. RXC J0616.8-4748 in the R band. Stars have been Figure A4. The centre of RXC J0006.0-3443 in the R band, 
excised from this image. Stars have been excised from this image. 


APPENDIX A: ADDITIONAL FIGURES 

Figure Figure |A2| Figure |A3| and Figure |A4| show R 
band images of a selection of our targets, where the stars 
have been excised from the images. 
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Figure Bl. Radial density profiles for RXC J0006.0-3443. The 
points with uncertainties represent the galaxy count density pro¬ 
file, normalised by r| 00 , and the red dashed line is the best fitting 
model with parameters given in Table [Bl] 



Figure B2. Radial density profiles for RXC J0049.4-2931. The 
lines are the same as described in Figure [bT] 

APPENDIX B: INDIVIDUAL CLUSTER 
RADIAL PROFILES 

Figures |B1| to |B14| show radial count density profiles and 
best fitting models, before background subtraction, of in¬ 
dividual clusters in the REXCESS sample. The best fit¬ 
ting radial profile parameters for each galaxy population in 
individual clusters in the REXCESS sample are given in 
Table IB 1 1 


Figure B3. Radial density profiles for RXC J0345.7-4112. The 
lines are the same as described in Figure [Bl] 



Figure B4. Radial density profiles for RXC J0547.6-3152. The 
lines are the same as described in Figure [Bl] 
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Figure B5. Radial density profiles for RXC J0605.8-3518. The 
lines are the same as described in Figure [Bl] 



Figure B6. Radial density profiles for RXC J0616.8-4748. The 
lines are the same as described in Figure [Bl] 



Figure B7. Radial density profiles for RXC J0645.4-5413. The 
lines are the same as described in Figure [Bl] 



Figure B8. Radial density profiles for RXC J0821.8+0112. The 
lines are the same as described in Figure [Bl] 
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Figure B9. Radial density profiles for RXC J2023.0-2056. The 
lines are the same as described in Figure [Bl] 



Figure BIO. Radial density profiles for RXC J2048.1-1750. The 
lines are the same as described in Figure [Bl] 



Figure Bll. Radial density profiles for RXC J2129.8-5048. The 
lines are the same as described in Figure [Bl] 



Figure B12. Radial density profiles for RXC J2218.6-3853. The 
lines are the same as described in Figure [Bl] 
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Figure B13. Radial density profiles for RXC J2234.5-3744. The 
lines are the same as described in Figure [Bl] 



Figure B14. Radial density profiles for RXC J2319.6-7313. The 
lines are the same as described in Figure [Bl] 


© 2015 RAS, MNRAS 000, 000-000 



















vi Holland et al. 

Table Bl. NFW fitting results. 


Object 

Galaxy filter 

p 

So 

r ~ 3 

'500 

SbgNFW 

r" 2 

7 500 

RXC J0006.0-3443 

Bright, red 

2.7 ± 0.5 

49.7 ± 9.0 

7.2 ± 2.1 

RXC J0049.4-2931 

Bright, red 

3.4 ±0.5 

59.0 ± 10.9 

1.9 ±0.5 

RXC J0345.7-4112 

Bright, red 

3.2 ± 1.0 

12.2 ±4.5 

2.3 ±0.6 

RXC J0547.6-3152 

Bright, red 

2.67 ± 0.35 

50.1 ± 6.8 

5.4 ± 1.4 

RXC J0605.8-3518 

Bright, red 

3.2 ±0.4 

55.5 ±8.8 

7.2 ±0.8 

RXC J0616.8-4748 

Bright, red 

2.58 ± 0.26 

30.1 ± 3.4 

2.7 ± 0.5 

RXC J0645.4-5413 

Bright, red 

2.74 ± 0.29 

96.2 ± 13.0 

9.8 ± 1.5 

RXC J0821.8+0112 

Bright, red 

2.5 ±0.4 

25.8 ±4.0 

1.8 ±0.9 

RXC J2023.0-2056 

Bright, red 

2.7 ± 0.7 

23.6 ±5.9 

0.2 ± 1.5 

RXC J2048.1-1750 

Bright, red 

2.3 ±0.5 

37.3 ±7.7 

5.8 ± 2.2 

RXC J2129.8-5048 

Bright, red 

3.3 ±0.7 

46.1 ± 10.7 

3.3 ± 1.4 

RXC J2218.6-3853 

Bright, red 

2.9 ± 1.0 

20.1 ± 6.4 

5.1 ± 1.1 

RXC J2234.5-3744 

Bright, red 

3.18 ± 0.15 

129.4 ± 6.9 

8.4 ± 0.8 

RXC J2319.6-7313 

Bright, red 

1.6 ±0.4 

4.5 ±0.9 

2.2 ± 0.5 

RXC J0006.0-3443 

Faint, red 

2.0 ±0.5 

20.5 ±3.6 

36.1 ± 3.4 

RXC J0049.4-2931 

Faint, red 

3.21 ± 0.30 

59.7 ± 7.0 

15.5 ±0.7 

RXC J0345.7-4112 

Faint, red 

2.8 ±0.9 

39.7 ± 12.5 

5.4 ± 3.1 

RXC J0547.6-3152 

Faint, red 

2.4 ±0.6 

33.5 ±7.6 

25.5 ±3.3 

RXC J0605.8-3518 

Faint, red 

2.8 ±0.6 

45.6 ± 10.3 

33.7 ± 2.1 

RXC J0616.8-4748 

Faint, red 

2.3 ±0.7 

22.3 ±6.6 

26.8 ±2.6 

RXC J0645.4-5413 

Faint, red 

1.6 ±0.5 

27.6 ±6.5 

44.0 ±4.6 

RXC J0821.8+0112 

Faint, red 

1.7 ± 0.5 

15.1 ± 3.2 

13.3 ±3.1 

RXC J2023.0-2056 

Faint, red 

2.1 ± 0.7 

16.7 ± 3.9 

12.1 ± 3.9 

RXC J2048.1-1750 

Faint, red 

2.3 ±0.5 

51.4 ± 10.0 

45.6 ±4.2 

RXC J2129.8-5048 

Faint, red 

2.8 ±0.5 

54.2 ± 8.9 

21.2 ±2.7 

RXC J2218.6-3853 

Faint, red 

3.9 ± 1.5 

36.0 ± 17.8 

29.9 ± 1.6 

RXC J2234.5-3744 

Faint, red 

2.6 ±0.7 

40.4 ± 10.2 

42.1 ± 4.0 

RXC J2319.6-7313 

Faint, red 

2.3 ± 1.3 

11.2 ±4.8 

17.8 ±2.1 

RXC J0006.0-3443 

Bright, blue 

1.7 ± 0.9 

4.9 ± 1.4 

5.9 ± 1.7 

RXC J0049.4-2931 

Bright, blue 

1.2 ± 1.2 

1.4 ± 0.9 

2.0 ±0.8 

RXC J0345.7-4112 

Bright, blue 

6.6 ±2.1 

(1.8 ± 1.0) x 10 2 

0.96 ±0.20 

RXC J0547.6-3152 

Bright, blue 

4.4 ±6.1 

(2.6 ±3.8) x 10 1 

6.3 ± 1.2 

RXC J0605.8-3518 

Bright, blue 

9.4 ± 13.0 

(1.4 ±2.5) x 10 2 

6.8 ± 0.5 

RXC J0616.8-4748 

Bright, blue 

2.9 ± 1.3 

7.3 ±3.5 

3.8 ±0.6 

RXC J0645.4-5413 

Bright, blue 

(0.1 ± 2.4) x 10 2 

(0.5 ±6.5) x 10 2 

9.7± 0.7 

RXC J0821.8+0112 

Bright, blue 

2.0 ±0.9 

5.2 ± 1.7 

0.8 ± 1.0 

RXC J2023.0-2056 

Bright, blue 

4.1 ± 2.5 

33.9 ±21.7 

0.4 ± 1.2 

RXC J2048.1-1750 

Bright, blue 

1.1 ± 0.5 

3.4 ± 1.0 

5.9 ± 1.1 

RXC J2129.8-5048 

Bright, blue 

7.5 ±3.3 

(1.6 ± 1.4) x 10 2 

2.47 ± 0.26 

RXC J2218.6-3853 

Bright, blue 

2.7 ± 0.9 

7.3 ±3.1 

3.6 ±0.5 

RXC J2234.5-3744 

Bright, blue 

1.0 ± 1.4 

2.4 ± 1.2 

11.9 ±3.7 

RXC J2319.6-7313 

Bright, blue 

5.3 ±2.6 

55.2 ± 34.1 

4.3 ± 0.5 

RXC J0006.0-3443 

Faint, blue 

-0.6 ±0.6 

0.99 ± 0.28 

38.2 ± 10.6 

RXC J0049.4-2931 

Faint, blue 

3.0 ±0.7 

20.4 ±5.8 

29.9 ±0.9 

RXC J0345.7-4112 

Faint, blue 

(2.0 ±5.0) x 10 1 

(0.0 ± 1.7) x 10 6 

7.7 ± 0.4 

RXC J0547.6-3152 

Faint, blue 

0.4 ±0.7 

2.6 ± 1.0 

44.5 ±4.7 

RXC J0605.8-3518 

Faint, blue 

1.8 ±3.2 

3.6 ± 5.1 

70.4 ± 4.4 

RXC J0616.8-4748 

Faint, blue 

-1.1 ± 1.1 

0.10 ±0.08 

37.3 ±2.7 

RXC J0645.4-5413 

Faint, blue 

0.6 ± 1.0 

1.5 ±0.8 

76.2 ±2.0 

RXC J0821.8+0112 

Faint, blue 

0.9 ±0.7 

3.6 ± 1.0 

7.9 ±3.6 

RXC J2023.0-2056 

Faint, blue 

4.4 ±3.9 

14.8 ±21.4 

6.9 ± 1.3 

RXC J2048.1-1750 

Faint, blue 

1.4 ±0.4 

11.7 ± 2.0 

56.9 ±2.6 

RXC J2129.8-5048 

Faint, blue 

-1.5 ±0.5 

0.078 ± 0.029 

0.7 ± 10.5 

RXC J2218.6-3853 

Faint, blue 

2.6 ± 13.3 

1.5 ±8.1 

47.6 ±3.2 

RXC J2234.5-3744 

Faint, blue 

4.4 ± 1.4 

70.2 ± 32.7 

108.4 ±2.5 

RXC J2319.6-7313 

Faint, blue 

1.5 ± 1.3 

4.0 ± 2.4 

33.3 ±2.9 
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APPENDIX C: BACKGROUND COUNT 
DENSITY ANALYSIS 

The background count densities for all of the objects, us¬ 
ing the four main galaxy population filters are shown in 
Table ICll 
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Table Cl. Background count density measurements. In the case of the RXC J2023.0-2056 bright blue filter, no objects are detected in 
the region used for measuring Sbgsimpie- 


Object 

Galaxy filter 

•^bg simple 
arcmin -2 

^Sbgsector 

arcmin -2 

‘S'bgNFW 

arcmin -2 

^bgsimple/'S'bgNFW 

RXC J0006.0-3443 

Bright, red 

0.166 + 0.020 

0.02 

0.12 + 0.04 

1.34 

RXC J0006.0-3443 

Faint, red 

0.70 ± 0.04 

0.04 

0.62 ± 0.06 

1.12 

RXC J0006.0-3443 

Bright, blue 

0.128 + 0.017 

0.02 

0.102 ± 0.029 

1.25 

RXC J0006.0-3443 

Faint, blue 

1.14 + 0.05 

0.07 

0.66 + 0.18 

1.73 

RXC J0049.4-2931 

Bright, red 

0.064 + 0.011 

0.01 

0.051 ± 0.013 

1.25 

RXC J0049.4-2931 

Faint, red 

0.427 ± 0.027 

0.03 

0.411 ± 0.019 

1.04 

RXC J0049.4-2931 

Bright, blue 

0.074 + 0.011 

0.01 

0.054 + 0.022 

1.38 

RXC J0049.4-2931 

Faint, blue 

0.81 ± 0.04 

0.04 

0.795 ± 0.024 

1.02 

RXC J0345.7-4112 

Bright, red 

0.029 + 0.010 

0.01 

0.027 + 0.007 

1.08 

RXC J0345.7-4112 

Faint, red 

0.095 ± 0.019 

0.02 

0.06 ± 0.04 

1.50 

RXC J0345.7-4112 

Bright, blue 

0.011 ± 0.006 

0.00 

0.0111 ± 0.0024 

0.98 

RXC J0345.7-4112 

Faint, blue 

0.098 + 0.019 

0.02 

0.089 ± 0.004 

1.11 

RXC J0547.6-3152 

Bright, red 

0.182 ± 0.024 

0.05 

0.134 + 0.034 

1.35 

RXC J0547.6-3152 

Faint, red 

0.71 ± 0.05 

0.06 

0.63 ± 0.08 

1.13 

RXC J0547.6-3152 

Bright, blue 

0.179 + 0.024 

0.05 

0.155 ± 0.030 

1.15 

RXC J0547.6-3152 

Faint, blue 

1.26 + 0.06 

0.22 

1.10 + 0.12 

1.14 

RXC J0605.8-3518 

Bright, red 

0.209 + 0.019 

0.03 

0.184 + 0.020 

1.13 

RXC J0605.8-3518 

Faint, red 

0.90 + 0.04 

0.06 

0.87 + 0.05 

1.03 

RXC J0605.8-3518 

Bright, blue 

0.180 + 0.017 

0.02 

0.175 ± 0.014 

1.03 

RXC J0605.8-3518 

Faint, blue 

1.84 + 0.06 

0.06 

1.81 ± 0.11 

1.01 

RXC J0616.8-4748 

Bright, red 

0.087 + 0.012 

0.01 

0.060 + 0.012 

1.44 

RXC J0616.8-4748 

Faint, red 

0.644 ± 0.033 

0.04 

0.60 + 0.06 

1.07 

RXC J0616.8-4748 

Bright, blue 

0.088 + 0.012 

0.03 

0.086 + 0.013 

1.03 

RXC J0616.8-4748 

Faint, blue 

0.99 + 0.04 

0.05 

0.84 ± 0.06 

1.18 

RXC J0645.4-5413 

Bright, red 

0.302 ± 0.022 

0.08 

0.231 ± 0.035 

1.31 

RXC J0645.4-5413 

Faint, red 

1.19 + 0.04 

0.07 

1.04 + 0.11 

1.14 

RXC J0645.4-5413 

Bright, blue 

0.225 ± 0.019 

0.08 

0.230 + 0.016 

0.98 

RXC J0645.4-5413 

Faint, blue 

1.86 + 0.06 

0.07 

1.81 ± 0.05 

1.03 

RXC J0821.8+0112 

Bright, red 

0.055 ± 0.011 

0.02 

0.031 ± 0.016 

1.77 

RXC J0821.8+0112 

Faint, red 

0.301 ± 0.027 

0.03 

0.23 + 0.05 

1.28 

RXC J0821.8+0112 

Bright, blue 

0.029 ± 0.008 

0.01 

0.013 + 0.017 

2.13 

RXC J0821.8+0112 

Faint, blue 

0.224 ± 0.023 

0.02 

0.14 + 0.06 

1.60 

RXC J2023.0-2056 

Bright, red 

0.011 ± 0.011 

0.01 

0.002 ± 0.014 

5.66 

RXC J2023.0-2056 

Faint, red 

0.12 + 0.04 

0.04 

0.106 + 0.034 

1.15 

RXC J2023.0-2056 

Bright, blue 

0.0 + 0 

0.00 

0.004 + 0.010 

0.00 

RXC J2023.0-2056 

Faint, blue 

0.066 + 0.027 

0.02 

0.061 ± 0.011 

1.09 

RXC J2048.1 1750 

Bright, red 

0.245 ± 0.022 

0.05 

0.16 + 0.06 

1.55 

RXC J2048.1 1750 

Faint, red 

1.37 + 0.05 

0.08 

1.24 + 0.11 

1.11 

RXC J2048.1 1750 

Bright, blue 

0.238 + 0.021 

0.03 

0.161 ± 0.029 

1.48 

RXC J2048.1 1750 

Faint, blue 

1.70 + 0.06 

0.06 

1.54 + 0.07 

1.10 

RXC J2129.8-5048 

Bright, red 

0.051 ± 0.014 

0.01 

0.038 + 0.017 

1.33 

RXC J2129.8-5048 

Faint, red 

0.290 + 0.034 

0.07 

0.248 + 0.031 

1.17 

RXC J2129.8-5048 

Bright, blue 

0.027 + 0.010 

0.01 

0.0289 ± 0.0031 

0.95 

RXC J2129.8-5048 

Faint, blue 

0.204 ± 0.028 

0.04 

0.01 ± 0.12 

24.40 

RXC J2218.6-3853 

Bright, red 

0.130 + 0.015 

0.02 

0.115 ± 0.025 

1.13 

RXC J2218.6-3853 

Faint, red 

0.673 + 0.034 

0.05 

0.678 + 0.035 

0.99 

RXC J2218.6-3853 

Bright, blue 

0.095 ± 0.013 

0.02 

0.081 ± 0.011 

1.18 

RXC J2218.6-3853 

Faint, blue 

1.11 ± 0.04 

0.04 

1.08 + 0.07 

1.03 

RXC J2234.5-3744 

Bright, red 

0.226 + 0.022 

0.02 

0.169 + 0.017 

1.34 

RXC J2234.5-3744 

Faint, red 

0.92 ± 0.04 

0.04 

0.84 ± 0.08 

1.10 

RXC J2234.5-3744 

Bright, blue 

0.288 + 0.024 

0.02 

0.24 + 0.07 

1.21 

RXC J2234.5-3744 

Faint, blue 

2.19 + 0.07 

0.06 

2.17 + 0.05 

1.01 

RXC J2319.6-7313 

Bright, red 

0.088 + 0.012 

0.01 

0.051 ± 0.010 

1.74 

RXC J2319.6-7313 

Faint, red 

0.431 ± 0.028 

0.02 

0.41 ± 0.05 

1.05 

RXC J2319.6-7313 

Bright, blue 

0.101 ± 0.013 

0.01 

0.100 + 0.011 

1.00 

RXC J2319.6-7313 

Faint, blue 

0.81 ± 0.04 

0.06 

0.77 + 0.07 

1.06 
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